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Abstract 

We survey some topics involving the Whitham equations, concentrating on 
the role of iptp* (or square eigenfunctions) in averaging and in producing Cauchy 
kernels and differentials on Riemann surfaces. 
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1 INTRODUCTION 



The Whitham equations and averaging methods have become a very important part 
of modern mathematical physics (see e.g. |l|, 
1431, ~ 



20), |25|, |26|, E8L |29|, |35|, |36|, |39|, |41 



52, 55, 



7C, 71, 72, 73, 82, 80, 33, 34, M 89, 90, m 



99], |100|| ). Most of the necessary background and details are to be found somewhere 



in these references but there are also a number of formulas and assertions which 
need clarification or amendment. We will look at some of this here and attempt to 
organize a subset of information in a coherent and more or less rigorous manner with 
the aim of "understanding" the basic ideas. Connections to dispersionless theory are 
indicated and it also seems desirable to indicate explicitly the main features that 
emerge which are relevant to work in topological field theory (TFT) etc., so we will 
display a lot of this. More recently there has also been a surge of activity connected to 
N = 2 supersymmetric gauge theories and some of this is sketched briefly as well (but 
certainly not explained). Thus the presentation was designed to be mainly expository 
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but the development also leads to some new connections and results. In particular 
we indicate various formulas for some Cauchy type kernels, give some connections 
of classical inverse scattering to averaging of square eigenfunctions for KdV, and 
provide some relations between averaged quantities for if>*if) and various differentials 
in general situations. In this spirit the classical Whitham theory can be phrased via 
ifjif)* which generates Cauchy type kernels as well and leads, in a sense indicated, 
to the dispersionless limit situation in genus zero. Relations between averaging and 
moduli spaces are developed following Dubrovin and Krichever and some Landau- 
Ginzburg (LG) type examples are displayed. Many basic formulas involving theta 
functions, differentials, etc. are also exhibited in order to make the text useful as 
a launching pad for further investigation. We would like to acknowledge valuable 
discussions with P. Grinevich, Y. Kodama, and A. Yaremchuk on related topics. 



2 SOME BACKGROUND 

We start from the viewpoint of Riemann surfaces. Thus (cf. [|7], || [28], |38|, |40], f|8 



58], |85| take an arbitrary Riemann surface E of genus g, pick a point Q and a local 
variable 1/k near Q such that k(Q) = oo. Let D = P\ + • • • + P g be a positive divisor 
of degree g and write ip for the (unique up to a constant multiplier) Baker- Akhiezer 
(BA) function characterized by the properties (A) if) is meromorphic on E except for 
Q where if)(P)exp(—q(k)) is analytic (q(k) = kx + k 2 y + k 3 t will do for illustration and 
(*) if) ~ exp(q(k))[l + X)i°(£j near Q) (B) On E/Q, if) has only a finite number 
of poles. In fact if) can be taken in the form (P £ E, Pq ^ Q) 

r f P , ~i ^2 Q(A(P) + xU + yV + tW + Zo) /f _ N 

J Pa <d(A{P) + z ) 

where fi 1 = dk + ■ ■ ■ , fi 2 = d(fc 2 ) + • ■ ■ , fi 3 = d(k 3 ) +---,Uj = f b . Q 1 , Vj = 
J b . Q 2 , Wj = J b . Q 3 (j = 1, • ■ ■ ,g), z = —A(D) — K, and a is a normalizing fac- 
tor (see below). The symbol ~ will be used generally to mean "corresponds to" or 
"is associated with" ; occasionally it also denotes asymptotic behavior and this should 
be clear from the context. Here the Qj are meromorphic differentials of second kind 
normalized via J a Qj = (aj, bj are canonical homology cycles) and we note that 
Q = xQ 1 + uQ 2 + tQ 3 ~ dq(k) normalized; A is the Abel-Jacobi map {A(P) = J P u)k 
where the u>k are normalized holomorphic differentials, k = l,---,g, J a .ujk = Sjk), 
and K = (Kj) ~ Riemann constants (2K = —A{Ky) where K a is the canonical 
class of E ~ equivalence class of meromorphic differentials) so Q(A(P) + z ) has ex- 
actly g zeros (or vanishes identically). The paths of integration are to be the same 
in computing Jp Q Q l or A(P) and it is shown in || |7j that if) is well defined (i.e. 
path independent). Then the £ 3 - in (*) can be computed formally and one determines 



3 



Lax operators L and A such that d y if> = Lip with dtip = Aip. Indeed, given the £ 3 - 
write u = —2d x ^i with w = 3£i<9 x .£i — 3d 2 £i — 3d x ^ 2 - Then formally, near Q, one has 
{-d y + d 2 + u)ip = 0(l/k)exp(q) and {-d t + d 3 x + {3/2)ud x + w)ip = 0(l/k)exp(q) 
(i.e. this choice of u, w makes the coefficients of k n exp(q) vanish for n = 0, 1, 2, 3 
note the similarity here to D-bar dressing techniques as in |1(J). Now define 



L = d 2 + u and A = d x + {3/2)ud x + w so d y if) = Lip and d t if) = Aip. This fol- 
lows from the uniqueness of BA functions with the same essential singularity and 
pole divisors. Then we have, via compatibility L t — A y = [A,L], a KP-1 equation 
(3/4)u yy = d x [u t — {1/ 4){Quu x + u xxx )] and therefore KP equations are parametrized 
by nonspecial divisors or equivalently by points in general position on the Jacobian 
variety J(S). For more in this spirit we refer to work on the Schottky problem and 
the Novikov conjecture in |9|, [MJ], [57], |93|. The flow variables x, y, t arise in (A) via 



q(k) and then miraculously reappear in the theta function via xU + yV + tW; thus 
the Riemann surface itself contributes to establish these as linear flow variables on 
the Jacobian and in a certain sense defines the flow variables. The pole positions Pj 
do not vary with x, y, t and (f ) u = 2d1logQ{xU + yV + tW + z ) + c exhibits 6 as a 
tau function (see |3C| for c). 



Now a divisor D* of degree g is dual to D (relative to Q) if D + D* is the null 
divisor of a meromorphic differential Q = dk + {(3/k 2 )dk + • • • with a double pole 
at Q (look at ( = 1/k to recognize the double pole). Thus D + D* — 2Q ~ 
so A{D*) - A{Q) + K = — [A(D) - A{Q) + K]. One can define then a function 

ip*(x,y,t,P) = exp(-kx - k 2 y - k 3 t)[l + £1/ k) H ] based on D* (dual BA function, 

7^ if)* - see Section 6.8), and a differential Q with zero divisor D + D*, such that 
= ipifj*Q is meromorphic, having for poles only a double pole at Q (the zeros of Q 
cancel the poles of if>if>*). Thus ipip*Q ~ ifiifi* (l + (/5/ k 2 + ■ • -)(ifc is meromorphic with a 
second order pole at oo, and no other poles. For L* = L and A* = — A + 2w — {3/2)u x 
one has then (d y + L*)ip* = and (c^ + A*)^* = 0. Now the prescription above seems 
to specify for ip* (U = xU + yV + tW, z* = -A(D*) - K) 

„ -J->n 1+ ,n 2+ tn 3 ) fl(.A(P) -U + z* ) 
V 0(A(P) + z* o ) [ - } 

Various further hypotheses on E lead to the traditional KP situations (KP-1 or 
KP-2) or to various reductions (e.g. nKdV). We mention again especially |25], 
|2"T| , P8] , |5B] , |5^ , |6^ , and also some recent work on 2-D periodic problems as 



in [^7J, [7J], [72]] . One is generally interested in real smooth solutions and in this 
respect KP-1 (~ Xij imaginary, x<i = y) presents problems - it is not even formally 
integrable in general (cf. |>7J for details). We will not discuss general theory 



here but always assume we have some integrable (finite zone) situation related to an 
appropriate Riemann surface (and perhaps periodic potentials). The analysis relative 
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to asymptotics and averaging will be cursory and formal at times, in keeping with 
"tradition" here, but we have tried to be somewhat more precise (there are of course 
some precise theorems in the literature but a discussion of these would often involve 
more hypotheses and detail than appropriate for this exposition). In particular we 
try to clarify some of Krichever's work which can be difficult to understand at times. 
Our interest is basically to sketch a framework, with formulas and results relative to 
Riemann surfaces, which is sufficiently detailed to be believable, and in which some 
meaningful subset of a full theory can be displayed. 



THE WHITHAM EQUATIONS AND AVER- 
AGING 



Averaging can be rather mysterious at first due to some hasty treatments and bad 
choices of notation - plus many inherent difficulties. Some of the clearest exposition 
seems to be in |], ||, |2(| [39], [41], |82|, |100|| whereas the more extensive developments in 

301 



e.g. |!5|, |26[ g|, |29|, |35|, |54|, |66|, |67 



89] , |90[| may be confusing at first, until one 



realizes what is going on. We will try to be careful and distinguish situations with 
periodic potentials from more general averaging situations involving ergodic ideas. 



3.1 Remark on PDE connections 

One choice of background situation involves an examination of dispersionless limits 
with no concern for periodicity or Riemann surfaces (cf. here || |15|, [H], |Tj|, 0). Thus 
given e.g. a KdV equation u t + 6uu x + u xxx = set ex = X and et = T leading to 
u T + 6u e u e x + t 2 u xxx = 0. The study of u e — > u where u satisfies the Euler equation 
Ut + Quu x = is a very delicate matter, of great interest in applications and in PDE 
but we do not discuss this here (cf. []79|). On the other hand the purely algebraic 
passage of the background mathematics of KdV (involving Lax operators, the KdV 
hierarchy, tau functions, vertex operators, etc.) to the corresponding background 
mathematics of the dispersionless theory, is relatively easy and will be indicated 
below. Moreover it is of great importance in an entirely different direction, namely 
in the study of topological field theory, strings, and 2-D gravity (see e.g. || O, [HI 



18, 24, ml B2L M m, 58, M, UR FA, PHI). We will insert such material later as 



appropriate. 



3.2 Preliminary ideas 



Let us follow H| (cf. also f67|) in order to have a suitably complicated example 
leading directly to matters of interest here, so as background we consider the KP 
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equation 

3 3 1 

-a 2 u yy = d x [u t - -uu x + -u xxx ] (3.1) 

(a 2 = 1 ~ KP-2; a 2 = -1 ~ KP-1). Here aL ~ fl£ - u with <70„V = ~ and 
A ~ d^ — (3/2)ud x + w, which is slightly different than before (to connect notations let 
u — ► — u and insert cr). Then one has the compatibility equations [d y —L, d t — A] = as 
before and it should be noted that for any function g(t), t ~ x,y, t, t 4 , • • •, operators 
L = gLg~ x + d y g g~ x and A = gAg~ x + d t g g~ x could be used, corresponding to a 
new wave function ip — gip. The notation of |58| also involves differentials Qj ~ Q l 



where Q { = dk\l + 0(l/k)} and J ah Qj = with Uf = (l/27ri) / 6fc (so Z7* ~ 
t/f ~ Vfc, and t/f ~ H4, up to factors of 2m). There will be solutions if) as before 



in ( |2.1| ) with potentials u given via (f). We recall that Riemann theta functions 
are more precisely written as Q(z\B) where B ~ a certain matrix (cf. below), and 
z ~ (21, • • • , z g ) ~ xf/ + + tW + zo for example. One knows via [§, |30|, [87], |94| that 
u of the form (f) u = 2d 2 logQ(z\B) + c is a solution of KP if and only if the matrix 
B defining is the b-period matrix of a Riemann surface determined via J b dujj 
(with z as defined). Now consider the spectral theory of 2-D periodic operators (**) 
(ad y — d 2 + u{x 1 y))ip = where u(x,y) = u(x + ai,y) = u(x, y + a 2 ). Bloch solutions 
are defined via 

ip(x + a t ,y, w u w 2 ) = wiij)(x, y, w b w 2 ); (3.2) 

ip(x, y + a 2 , u»i, w 2 ) = w 2 ^Jj(x, y, w 2 ) 

and we assume ip(0,0,wi,w 2 ) = 1. The pairs Q = (wi,w 2 ) for which there exists 
such solutions is called the Floquet set T and the multivalued functions p(Q) and 
E(Q) such that w% = exp{ipa\) with w 2 = exp{iEa 2 ) are called quasi-momentum 
and quasi-energy respectively. By a gauge transformation ip — > exp(h(y))ip, with 
d y h(y) periodic one obtains solutions of (**) with a new potential u = u — ad y h so 
we can assume /q 1 u(x, y)dx = 0. 



3.3 Floquet theory 

For M = ady — d 2 (with u = 0) the Floquet set is parametrized by k G C such that 
= expiikai) and w° — exp(—k 2 a 2 /a) and the Bloch solutions are i[)(x,y,k) = 
exp(ikx — k 2 y/a). The adjoint Bloch solutions are ip + (x, y, k) = exp(—ikx + k 2 y/a) 
satisfying (ad y + d 2 )ip + = 0. The image of the map k — ► w°) € C 2 is the Floquet 
set for Mo corresponding to the Riemann surface with intersections corresponding to 
pairs k 7^ k! such that wf(k) = wf(k'), i — 1,2. This means k — k! = (2nN/ai) 
and A; 2 - (A;') 2 = (2maM/a 2 ) where iV, M are integers, creating "resonant" points 
k = kw,M = (jN/ai) — (icrMai/Na 2 ), N 7^ 0, k! = k^N-M- Then for k 7^ fcw.M 
and it sufficiently small one can construct a formal Bloch solution of (**) in the form 
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of a convergent perturbation series (for any a). Thus outside of some neighborhoods 
of the resonant points one can obtain a Bloch solution ip(x, y, k Q ) which is analytic 
in ho, but the extension of ip to a resonant domain can be tricky. For 3?cr = (as 
in KP-1) the resonant points are dense on the real axis whereas for dta 7^ (as in 
KP-2) there are only a finite number of resonant points in any finite domain of C. 
In the latter case one can glue handles between domains around resonant points and 
create a Riemann surface T of Bloch solutions ip(x,y,Q), Q G T. Moreover, if the 
potential u can be analytically extended into a domain \Qx\ < T\, \^sy\ < t 2 , then 
the technique can also be adapted even when u is not small. Such Bloch solutions, 
normalized by ip(0, 0, Q) = 1, are meromorphic on F and in the case of a finite 
number of handles a one point compactification of F is obtained so that ip is in fact 



the BA function for F. Generally speaking finite zone situations as in ( |2.1| ) with 
potentials given by Riemann theta functions are quasi-periodic in nature. To single 
out conditions for periodicity one asks for meromorphic differentials dp and dE on F 
having their only singularities at Q ~ point at 00 of the form dp = dk(l + 0(k~ 2 )) 
and dE = ia^ 1 dk 2 (l + 0(k" z )), normalized so that all periods are real, and satisfying, 
for any cycle C on F, § c dp = (27in c /ai) with § c dE = (27rm c /a 2 ) where ric, m c are 
integers. Then the corresponding potentials u(x,y) will have periods ai and a 2 in x 
and y, with multipliers W\(P) = exp[ia\ J P dp) and w 2 (P) = exp(ia,2 J P dE). We go 
now to finite zone (or quasiperiodic) situations as in (^j]) with potentials as in (f), 



where z ~ xU+yV+tW+z can be written as (z k ) = [x f b Vt}+y J b £l 2 +t f b fi 3 +z°] = 
(Cfc + z k)i k = 1) ' ' ' i9- Since 6(z + 2niN) = Q(z) we could set ( = i6 so that as 
a function of 9 = (8k), u is periodic of period 2ir in each variable 8k- Now one 
wants to consider modulated finite zone situations where solutions are of the form 
u = u (xU +yV +tW\I) = u (8i, ■ • ■ , 9 g \Ii, ■ ■ ■ ,I n ) where u is periodic in the 6j with 
U,V,W = U,V,W(I). One assumes there will be slow variables X = ex, Y = ey, 
and T = et with fast variables x,y,t so that "asymptotic" solutions 

u = u (-S(X, Y, T)\I(X, Y, T)) + e Ul (x, y, t) + e 2 u 2 (x, y,t) + --- (3.3) 
e 

can be envisioned. In practice the parameters Ik depend on the moduli of our Riemann 
surface, which are then allowed to change smoothly with the slow variables X, Y, T (so 
U, V, W also can change), and one looks for solutions (|3.3j) with uniformly bounded 



U\. We do not look for convergence of the series in (3^) nor check any other features 
of "asymptotic solution" or "asymptotic series". Such procedures are standard in the 
study of what are called weakly deformed soliton lattices. 
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3.4 KdV averaging 



We will motivate this study with examples from KdV, where one will see explicitly the 
nature of things. First from |82| , in a slightly different notation, write q t = 6qq x — q xxx 
with Lax pair L = —d 2 + q, B = — 4<9ijl + 3(qd x + d x q), L t = [B, L\, Lip = Xip, and 
ipt = Bip. Let ip and be two solutions of the Lax pair equations and set \I/ = ipcp; 
these are the very important "square eigenfunctions" which arise in many ways with 
interesting and varied meanings. Evidently \l/ satisfies 

[-dl + 2(qd x + d x q)]V = 4\d x V; d& = -2q x ^> + 2(q + 2\)d x V (3.4) 

From ( |3.4| ) one finds immediately the conservation law (C): d t [^} + d x [6(q — 2A)\I/ — 
2dm = 0. If one looks for solutions of of the form *(ar, t, A) = l+T,?[&j(x, t)]\- j 
as A — >• oo then one obtains a recursion relation for polynomial densities 

d x * j+ i = \-\dl + (qd x + (j = 1, 2, • • •); *o = 1 (3-5) 

Now consider the operator L = —d 2 x + q in L 2 (— oo, oo) with spectrum consisting of 
closed intervals separated by exactly N gaps in the spectrum. The 2N + 1 endpoints 
A/; of these spectral bands are denoted by — oo < Ao < Ai < ■ • • < \2N < oo and 
called the simple spectrum of L. They can be viewed as constants of motion for KdV 
when L has this form. We are dealing here with the hyperelliptic Riemann surface 
determined via -R 2 (A) = llo — ^k) an d one can think of a manifold Ai of iV-phase 
waves with fixed simple spectrum as an iV-torus based on 6j e [0, 2ir). Hamiltonians 
in the KdV hierarchy generate flows on this torus and one writes q = qjsi{B\, • • • , 9^) 
(the 9 variables originate as in our previous discussion if we use theta functions for the 
integration - cf. also below). Now there is no y variable so let us write 9j = xkj +tWj 
(we will continue to use du)j for normalized holomorphic differentials) For details 
concerning the Riemann surface we refer to |3|, ^] and will summarize here as follows. 
For any qjy as indicated one can find functions \ij[x y t) via t, A) = Yli (A— t)) 
where fij(x,t) G [\2j-1, ^2j] and satisfies 

d x [Mj = -2z(i?(/i j )/ JJ(/Xj- - At*)); (3.6) 

i¥=3 

2N 

d&j = -2i[2(E X k -2jft)] • (i?0*i)/ nO*i - to) 

ijtj i^j 

In fact the fij live on the Riemann surface of R(X) in the spectral gaps and as x 
increases fij travels from \2j-1 to \%j on one sheet and then returns to \2j-1 011 
the other sheet; this path will be called the j //-cycle (~ aj). In the present con- 
text we will write the theta function used for integration purposes as 0(z,t) = 
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J2 m £z N exp[iii (2 (m,z) + (m, mi))] where z £ Z N and r denotes the N x N period 
matrix (r is symmetric with 3r > 0). We take canonical cuts a{, b{ (i — 1, • - • , iV) 
where % ~ (\2j-2, ^2j-i) as in |3^] (where a picture is drawn). Let dujj be holo- 
morphic diffentials normalized via J a . duj^ = 5jk (the cycle cij corresponds to a loop 
around the cut a,). Then q^ can be represented in the form 

2N 

q N (x, t) = A + T - 2d?logQ(z(x, t); r); A = £ A 5 -; (3.7) 



7" = (7y) = (* dwj-); 7^ = -r y ; r = -2^ f Ac% 

Jbi 1 a J 

and z(x,t) = — 2i[c N (x — x ) + 2(Ac Af + 2c N ~~ 1 )t] + d where (c ff )i = qjv arises from 
the representation da;, = CijA- 7_1 )[rfA/i2(A)] (d is a constant whose value is not 
important here). Then the wave number and frequency vectors can be defined via 
k = —AiitT~ l c N and w = -8i7rr~ 1 [Ac Ar + 2c Ar - 1 ] with 9j(x,t) = KjX + Wjt + 6® (where 
the 9°j represent initial phases). 

To model the modulated wave now one writes now q = q^i^i, ■ ■ • , On] A) where 
Xj ~ Xj(X,T) and A ~ (A-,). Then consider the first 2A^ + 1 polynomial conservation 
laws arising from (|3.4|) - (|3.5|) and C for example (cf. below for KP) and write these as 
d t Tj(q) +d x Xj(q) = (explicit formulas are given in roughly as follows). We note 
that the adjoint linear KdV equation (governing the evolution of conserved densities) 
is d t jj + d^jj - 6qd x jj = (7,- ~ VfT,-) and Q3T5| ) has the form dj j+ i = (-(l/2)<9 3 + 
qd + dq) r yj. One then rewrites this to show that Qqd x ^j = d x [Qjj + i — Qq-^j + 3<9 2 7j] so 
that the adjoint equation becomes 

d tlj + d[-2d 2 l3 + Q qij - 6 7j+ i] = (3.8) 

which leads to fl3T9|) and fl3.13|) below (after simplification of Q3.8|) ). Then comes the 
crucial averaging step. Keep the slow variables X, T constant and average over the 
fast variable x to obtain 



d T <T j (q N )>+d x <X j (q N )>=0 (3.9) 
(note e.g. d t = edj>). The procedure involves averages 

1 f L 

< Tj(q N ) >= UrriL^oo— J ^Tj(q N )dx (3.10) 

for example (with a similar expression for < Xj(q^) >) and an argument based on 
ergodicity is used. Thus if the wave numbers Kj are incommensurate the trajectory 
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{qisr(x,t)] x G (—00,00)} will densely cover the torus M.. Hence we can replace x 
averages with 

1 rl-K r2ir _ N 

<^(^)>=7^yv/ r j { qN {e))j{de j (3.11) 



1 



For computational purposes one can change the 9 integrals to \x integrals and obtain 
simpler calculations. By this procedure one obtains a system of 2N + 1 first order 
partial differential equations for the 2N + 1 points Xj(X,T), or equivalently for the 
physical characteristics (k(X,T),w(X,T)) (plus < q^ >)• 

The above argument may or may not have sounded convincing but it was in 
any case very loose. Let us be more precise following |j3~9| . One looks at the KdV 
Hamiltonians beginning with H = H(q) = iim i _ DO (l/2.L) J L L (q 2 + (l/2)q 2 x )dx (this 
form is appropriate for quasi-periodic situations). Then q t = {q,H} where {/, g} = 
lim L ^ 00 (l/2L) f L L (5f/5q)d x (5g/5q)dx (averaged Gardner bracket). The other Hamil- 
tonians are found via 

d 8H^ = _ l^Hrn (m > Q); SHo = x 

oq I oq oq 

where jj ~ Vi?j ~ (5Hj/Sq) (cf. here || |39|). It is a general situation in the study 
of symmetries and conserved gradients (cf. UTI] ) that symmetries will satisfy the 
linearized KdV equation (d t — 6d x q + d x )Q = and conserved gradients will satisfy 
the adjoint linearized KdV equation (pt — Qqd x + d x )Q^ = 0; the important thing to 
notice here is that one is linearizing about a solution q of KdV. Thus in our averaging 
processes the function q, presumed known, is inserted in the integrals. This leads 
then to 

m = m = -231^ - 6^1 + 6,^ (3.13) 
with (|3.9|) holding, where < d 2 <p >= implies 

< X j >= lim^^r f L (-6^±i + 6q N 5 -^)dx (3.14) 
2L J-l oqN oqN 

Note that d x appears in Xj f° r example but in ( |3.9|) d x — > edx- This is somewhat 
confusing but apparently gives correct first order terms in e, after which e is cancelled 
out and eventually allowed to approach 0. In 3J| the integrals are then simplified in 



terms of fx integrals and expressed in terms of abelian differentials. This is a beautiful 
and important procedure linking the averaging process to the Riemann surface and 
is summarized in [[82] as follows (see also below for more details and a KP version). 
One defines differentials 

1 N rl\ 
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where the Cj, Ej are determined via § b . Oi = = & 2 (i = 1, 2, ■ • • , N). Then it can 
be shown that 

<*>~<r>~i:^;<*>~i:^ ("6) 

with ~< T > (^/e 2 ) and < X > (d£/£ 2 ) ~ 12[(d£/£ 4 ) - where // = -> 

oo (// ~ (l/v^) 1/2 ) so d// = -2£- 3 d£ {d£/£ 2 ) (£/2)dfi (dfi/2y/Ji). Since 

Cl t = 0{fi N /fi N+ ^)dfi = Ofa-Wdfi, Q 2 = 0{^ 2 )dfi (with lead term -(1/2)) we 
obtain < * >~< T >= 0(1) and < X >= 0(1). Thus (0), Q, (C) generate all 
conservation laws simultaneously with < Tj > (resp. < Xj >) giving rise to Oi (resp. 
Cl 2 ). It is then proved that all of the modulational equations are determined via the 
equation 

«9 t j\ = I2d x & 2 (3.17) 

where the Riemann surface is thought of as depending on X, T through the points 
Xj(X, T). In particular if the first 2iV + l averaged conservation laws are satisfied then 
so are all higher averaged conservation laws. These equations can also be written di- 
rectly in terms of the Xj as Riemann invariants via dxXj = SjdxXj for j = 0, 1, ■ • • , 2N 
where Sj is a computable characteristic speed. Thus we have displayed the prototyp- 
ical model for the Whitham or modulational equations. 

3.5 Extension to KP 



Given some knowledge of symmetries as sketched in [11] for example one is tempted 



to rush now to an immediate attempt at generalizing the preceeding results to finite 
zone KP via the following facts. The KP flows can be written as d n u = K n (u) where 
the K n are symmetries satisfying (in the notation of |pT| ) the linearized KP equation 
d 3 f3 = (l/A)d 3 (3+3d(u(3) + (3/A)d- 1 d 2 (3 = K'[(3\. The conserved densities or gradients 
7 satisfy the adjoint linearized KP equation = (l/4)<9 3 7 + ?>ud^ + (3/4)<9 -1 <9 2 7. 
Then, replacing the square eigenf unctions by ipip* one has e.g. ipip* = X)o° s nX~ n 
where s n ~ a 7„. Further d n u = K n+ i = ds n+ i = dRes L n = dVI n +i where V/ ~ 
Sf/Su (Res L n = nH*_i is generally used in the multipotential theory). We are 
working here in a single potential theory where all potentials «j in L = d+J2T u i+i^ % 
are expressed in terms of u 2 = u via operators with d and d^ 1 . One uses here 
the Poisson bracket {f,g} = J J(Sf/Su)d(Sg/Su)dxdy (Gardner bracket). Let us 
retrace the argument from [^] and see what applies for KP. Thus one has s n+ i ~ 
7 n+ i ~ V-Z„ + i as conserved gradients satisfying the adjoint linear KP equation (**) 



11 



<9t7 = (l/4)<9 3 7 + 3ui97 + (3/4)d~ 1 dy'j. The nonlocal term <9 _1 here could conceivably 
change some of the analysis. We need first a substitute for (|3.5|) or else a direct way 
of rewriting the adjoint equation as dt[A] + d x [B] = or perhaps we want rather 
d t [A] + d x [B] + d y [C] = 0. To get such a formula we can simply differentiate the 
adjoint KP equation to get 

dt{dl] + <9[i<9 3 7 + 3u<9 7 ] + y [~0y7] = (3.18) 

This removes all of the nonlocal terms and one doesn't have to deal with d[(3 / '4)<9 _1 dfcj] 

for example. Thus imagine a finite zone situation and write (7^ ~ V/j) 

1 . 3 
D,- = <9 7j ; ^ = -9 3 7j + 3ud^; Q 3 = -d ylj (3.19) 

so (|3.18fj ~ d t [Dj] + d[!Fj] + d y [Qj] = 0. One linearizes around a fixed finite zone 
solution u in the adjoint linear KP equation and puts this u into the Ij etc. Then 
average in (|3.18|) over the 9i variables as before to obtain 

d T < Vj > +d x < Tj > +dy < Gj >= (3.20) 

Then modeled on KdV one expects the terms in ( p.20[ ) to be expressible in terms 
of differentials but there is remarkably little information in this direction. First we 
expect < <9 7 j >= and < d y jj >= if e.g. d y ^j = ^Y^djj / d9i) and d^j = 
J2Xi(djj/ 'd6i) (see the analysis below). Then dx < >= ensues which is not 
very thrilling (and this would be the situation indicated by (|3.25|) - ( [3.26D below). 
Let us try then 

dth\ + di-^ + 30- x «fry] + dy^dyd-H = (3.21) 

Here we expect < 9 2 7 >= =< dyd^ 1 ^ > so set £ = ^d^iud^) with Cj = 
Sd^^udjj) and deduce that 

d T < 7 > +d x < £ >= 0; d T < 7,- > +d x < Cj >= (3.22) 

Then one wants to express < 7 > (resp. < jj >) and < C > (resp. < Cj >) 
in terms of differentials (and since < 7 >=< > one has some connection to 

Section 3.4). However order of magnitude considerations (see below) suggest that 
< <9 2 7 >= =< d y d~ 1 ^ > must be false - due to growth or whatever - and we take 

< t >=< i<9 2 7 + 3<T 1 («od7) > (3-23) 
with Q = (3/A)d y d- 1 j. Then 

d T < 7 > +d x < C> +d Y < Q >= (3.24) 

which conceivably might be useful. These matters will be covered after we develop 
another approach. 
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3.6 Averaging with ip*ip a la Krichever 

Let us therefore look at but in the spirit of (cf. also [0) which is more 
carefully done at times. We will expand upon this with some modifications in order 
to obtain a visibly rigorous procedure. Thus consider KP in the form Q3.1J ): 3a 2 u yy + 
d x (4u t — 6uu x + u xxx ) = via compatibility [d y — L,d t — A] = where L = <r _1 (<9 2 — u) 
and A = d 3 — (3/2)w<9 + u> (cr 2 = 1 is used in [[42] which we follow for convenience but 
the procedure should work in general with minor modifications - note d means d x ). We 
have then (d y — L)ip = with (pt — A)ip = and for the adjoint or dual wave function 
ip* one writes in |66j ip*L = —d y tp* with tp*A = d t ip* where ip*(fd j ) = (—dy(ip*f). 
The explicit formulas used in (and [J42[| ) are 

^ = eP x+Ey+nt+YX a nit . ^ Ux + Vy + Wt + C, P) (3.25) 

^ = ^ps-Ey-ftt-Y* 1 *ib . -Vy-Wt-t, P) (3.26) 

where Oi = cXj(P), p = p(P), E = E(P), Q = tt(P), etc. (see below for more detail). 
We will assume these to be correct although this is not really obvious since the proof 
in |66] is rather "heuristic" and the form is not immediately compatible with (|2.1| ). 
However the arguments to follow are essentially independent of this choice of notation. 
Now one sees immediately that 

= i>*14 + d x {i)*L l i)) + d 2 x {^L 2 ^) + ■■■ (3.27) 

where e.g. U = (—l) r /r\)(d r L/d(d) r . In particular 

r xx ^ = w>** - 2d(w x ) + d 2 (r*P) (3.28) 

- V xxx ^ = W*** - d\r^) + 3[<9 2 (V*^) - <WU] (3.29) 
This means (L* = L, vp*A ~ A*ip* = -d 3 i)* + (3/2)d(ijj*u) + wip*) 

(iJ)*L)iJ) = [(d 2 - u)ij)*\il) = i)* x J) - uip*ip = ip*Lip - 2d(ip*ip x ) + d 2 (ip*tlj) (3.30) 

which implies L 1 = —2d and L 2 = 1. Next 

(rAw = -r xxx ^ + + ru x )*p + = r w) + (3.31) 

+d[ip*(-u - 3d 2 )^j] + 3<9 2 (</>*^) - d 3 (tp*ij) 

so that A 1 = —3d 2 + (3/2)u, A 2 = 3d, and A 3 = —1. We think of a general Riemann 
surface S s . Here one picks holomorphic differentials duj^ as before and quasi-momenta, 
quasi-energies, etc. via dp ~ Q\, dE ~ Q 2 , dQ ~ Q3, ■ ■ ■ where A ~ k, p = Jp o Q\, 
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etc.), and use ~ C/| ~ and f/f ~ (general ~ Q k are discussed at 
various places in the text). Normalize the flk now so that 9ft J a . Qk = = 9ft J* 6 . fife 
and take then e.g. C/ ~ , — (l/27rz) J am (A;, m = 1, • • • , so that U, V, W are 
real 2g period vectors (this is essentially equivalent to previous normalizations, e.g. 
J a Qk = 0). Then one has BA functions ip(x, y, t, P) as in (ET25D or (gjg) . As before 
we look for approximations as in ( [3.3|) based on uo(xU + ?/V + iW|7) = u (9j,Ik)- 
For averaging 0j ~ xUj + yVj + tWj + Q, 1 < j < 2g, with period 2n in the 9j 
seems natural (but note 9j, 9 g+ j ~ Uj, etc. - cf. ( p.99| )). Then again by ergodicity 



< > x = lim L _ 00 (l/2L) f^ L (fidx becomes < >= (l/(2n) 2g J---J <j)d 2g 9 and one 
notes that < d x (p >= automatically for bounded. In jEJ one thinks of <j)(xU + ■ ■ •) 
with <fi x = E U i (d(j)/d9 i ) and / • • • J(d(j)/d9 i )d 2g 9 = 0. 

Now for averaging we think of uq ~ uq(-S\I) as in ( |3.3|) with S,I ~ S, I(X, Y, T), 
dx*? = Z7, cVS* = V, and c^S 1 = W. We think of expanding about u with d x — ► 
(9 X + e<9x- This step will cover both x and X dependence for subsequent averaging. 
Then look at the compatibility condition (f ) : dtL — d y A + [L, A] = 0. As before 
we will want the term of first order in e upon writing e.g. L = L Q + eLx + • ■ ■ and 
A = A + eAi + • ■ ■ where L , A are to depend on the slow variables X, Y, T. As 



indicated in |)6| the term [L, A] — > {L, A} where {L, A} arises upon replacing d x be 
d x + edx in all the differential expressions and taking the terms of first order in e. 
However the formula in [p6] is unclear so we compute some factors explicitly. In fact, 



according to |42| , one can write now, to make the coefficient of e vanish 

d t Lx-d y Ax + [L ,Ax} + [Li,A )+F = 0; F = drL-dyA+^dxA- A l d x L) (3.32) 

Thus F is the first order term involving derivatives in the slow variables. To clarify this 
let us write ~ d/dX) L e = (d+ed) 2 - (u + eux + - ■ •) = d 2 -u + e(2dd- Ul )+O(e 2 ) 
and A e = (d + ed) 3 - (3/2) (u + eu x + •■■)■ (d + ed) + w + ew x + ■■■ = d 3 - (3/2)u d + 
w + e(3d 2 d-(3/2)uxd-(3/2)u d+wx) + O(e 2 ). Write then A e = A + Axe+O(e 2 ) and 
L e = L + Lxe + O(e 2 ) with Ax = 3d 2 d- {3/2) Ul d- (3/2)u d + wx and Li = 2dd-u x . 
Note that L 1 = -2d ^ t\ and A 1 = -3d 2 + (3/2) w ^ Ax but we can write 

A x = -A l d - {3/2) Ul d + wx = -A 1 ^ + A x \ (3.33) 

Lx = -l}d-u x = -L l d + Lx 
with Lx and Ax as in [fSq] . Then (f) becomes 

d t L -d y A^L ,A }+t{d t Lx-d y Ax + d T L-dYA^L ,Ax]HLi,A,]}+O{t 2 ) (3.34) 

But for uq, w functions of the slow variables only, the d t L , d y A , [L ,A ] terms 
vanish and we note that 

[L Q , Ax] + [Lx, Aq] = [L Q , -A l d] + \—L l d, A ] + [L , A x ] + [Lx, A ] = (3.35) 
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= [L , At) + [Li, Aq] + A l dL - tfdAo - L Q A l d + A l}d 

Then one notes that dxL Q ~ dxL and <9x^4o ~ <9x^4 (dx ~ 9) so dropping the terms 
in (|3.35 ) with an inoperative d on the right we obtain (|3.32|) . Next one writes, using 



dtWL^) - d y (ip*Aiip) = r{L u - A ly + [L x , A] + [L, A^ifj 



(3.36) 



= ip*(d t Li - d y A x + [Lo, A,] + [L x , A ))if; + d x {- ■ ■) 

and via ergodicity in x, y, or t flows, averaging of derivatives in x, y, or t gives zero so 
from ( |3.32| ) and ( |3.36| ) we obtain the Whitham equations in the form < ip*Fip >= 
(this represents the first order term in e - the slow variables are present in L , A Q , if), 
and if)*). In order to spell this out in one imagines X, Y, T as a parameter £ and 
considers etc. (in their perturbed form) with 



o p(S)i+s(0y+n($)t+<T-C«) 



<j>(U(£)x + V(£)y + W{i)t + C(tmO) (3-37) 



where J2i g &id ~ cr • C an d ^* = exp(-px — Ey — fit — a ■ ()<fi*(—Ux — Vy - VF£ - Cl-0 
(no £ variation - i.e. assume p, E, fl, U, V, W, I fixed). We recall that one expects 
Afc = Afc(X, y, T) etc. so the Riemann surface varies with £. Also recall that x,y,t 
and X, Y, T can be considered as independent variables. Now as above, using (|3.27| ), 
we can write 



Note also from (|3T37D for P and C(0 fixed (0 ~ xU + + yW + C) 

%/>*V>(Ole=o = (pz + % + + 



(3.38) 



(3.39) 



where / ~ df/d£. In ||2| one assumes (without discussion) that it is also permitted 
to vary £ and hold e.g. the Ik constant while allowing say the P to vary. Now 
differentiate the left side dt(if>*ip(£)) of (|3.38| ) in £ and use (|3.39| ) to obtain 

%[<WW))]le=o = W*i/> + W ■ + ( 3 - 4 °) 
+{{px + Ey + tlt)d t {if)*ip) + (Ux + Vy + Wt) ■ d t {(f)*(f) e ) + / • d^fa)} 
Fixing x,y,t and averaging ( |3.4C ) in the 9 variables yields 

< dz[d t (r^(0)}k=a >= « < « >+W< > (3.41) 
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Next one differentiates the right side of ( |3.38| ) in £, averages, and equates to ( |3.41| ) 
to obtain 

(3.41) =< rM >-< 9 e [4(f^V(0)]lf=o >= (3.42) 
< ip*Aijj > -p < %l)*A l il) > -U- < il)*A l ij)e > 

We note here as in ( |3.39| ) - ( |3.4U| ) (the result is unaltered if one disregards £ dependence 
in A 1 , as well as A 1 action on the isolated x term, and assumes I fixed, since only 
more terms d x (- ■ •) would occur whose average vanishes) 



(3.43) 



(3.44) 



9«9 € (V*AV)k=o = d x [ip*A 1 (px + Ey + fii)V> + 

+iP*A 1 (Ux + Vy + Wt) ■ ij} 9 \ = P4)*A 1 ^ + U ■ ijfA 1 ^ 
+ (px + Ey + (ly)d x {%l)*A l i)) + (fix + Vy + Wt) ■ d x {ij)* A 1 ^) 
and ( 3.42|) follows. Rewriting (|3.42|) with £ ~ Y we obtain (note ip*ipe = <fi 

- < ^dyAtp >= -d Y tt < ip*ip > -d Y W < <j)*(j)e > - 

-dyp < f AV > -dyU- < ^Aty > 
Similarly one uses (cf. (|3.38| )) 

d y (rm) = - lmo - wlV(O) + • o (3.45) 

with £ ~ T to get 

< ip*d T Lip >= d T E < > +d T V- < <p*<p e > + (3.46) 
+d T p < V»*^V > +d T U- < V>*^Vfl > 

Finally note that 

dyirA'm) - dtWLWt)) = (3.47) 

= r[L\A(o -a)- A\m - Dm) + d x {. . •) 

Using this with £ ~ X one gets then as above 

< ^(tfdxA - A l d x L)^ >= d x Vt < V*^V > ~d x E < ^A 1 ^ > + 

+ d x W- < ^L 1 ^ > -d x V- < ip*A l tp e > (3.48) 
We recall that dxS = U, dyS = V, and dxS = W so there are compatibility relations 

d Y U = d x V; d T U = d x W; d T V = d Y W (3.49) 
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Now add up ( ggg ) and (|3^6| ), and subtract ( |OS| ) to get 

< V*[^ T L - + L l d x A - A 1 d x L]i) >=< ip*F^ >= = (3.50) 

= (d T E - d Y ty < > +(d T p - d x ty < > +(<9 X P - d Y p) < ^A 1 ^ > 

We observe that if one lets the point P on £ vary with £, while holding and in- 
fixed, then 

9f(^XO)lf=o = + J/d-B + tdtt)<f)*(f) (3.51) 
which, together with ( |3.38| ) and ( |3.45|) , yields e.g. 

dfi < 0*0 >= -dp < ip*A l ip >■ dE < 4>*4> >= -dp < ^*LV > (3.52) 

(cf. (|3~42] ) - O)). It follows then from (ET5CD and (gj|) that 

= (fi y - P T )dp + (pr - fi x )dP + (E x - py)rffi (3.53) 

Thus one arrives at a version of the Whitham equations in the form 

p T = Q x ; p Y = E x ; E T = Q Y (3.54) 

where the last equation establishes compatibility of the first two via pty — Pyt — 
Hjy — Ext — dx{^iy — Er) — 0. We feel that this derivation from [42| is important 
since it again exhibits again the role of square eigenfunctions (now in the form ijj*ifj) 
in dealing with averaging processes. In view of the geometrical nature of such square 
eigenfunctions (cf. |0, [13] for example) one might look for underlying geometrical 
objects related to the results of averaging (cf. here |I4[). Another (new) direction 
we will discuss later involves the Cauchy kernels expressed via and their disper- 



sionless limits (cf. |I5], [49], ^(], |5l], |91|). Connections to Section 3.5 are also discussed 
later. 



3.7 Relations to moduli spaces 



Now following |68], |74] let M. 9 n be the moduli space of smooth algebraic curves 
T g of genus g with local coordinates k~ l (P) in the neighborhoods of iV punctures 
P a (^o 1 (-Po) = 0) so M. g N ~ P a , k~ l , a = 1,---,N} (see Section 7.3 for a 
discussion in the context of LG theory and Hurwitz spaces of the moduli space M 9 n ^ 
M. 9 n). First look at g — (r ~ P 1 ) and for simplicity take N — 1, (Pi ~ oo) with 
k x (P) = P + J2Tv s P~ s , so M ,i ~ {v s , s = l, - }. Set A = {A { = (l,i)} and define 
meromorphic functions 



n<(p) = fi M (P) = J2 w hi,s ps = k\(P) + oiK 1 : 



(3.55) 
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(P ~ point or local coordinate here in a flagrant abuse of notation - see below). 
If we keep g = but allow more punctures with k a (P) = J2°?i v a,s(P — Pa) s then 
one defines A = {A = (a,i), a = 1, • • • , N, i = 1, 2, • • • ; for i = 0, a ^ 1} and 
■Mo,n = {P a ,Va,s} with meromorphic functions (O^j as in ( |3.55|) ) 

tta,i( p ) =J2 w ^s(P - P*)~ s = K(P) + 0(1); a*,i(oo) = (a ^ 1) (3.56) 

-i 

One can write then in fact 

n a AP) = ^-.l K ^ )dZa ; n a>0 (P) = -log{P - P a ) (a ? 1) (3.57) 

where C a is a small cycle around P a . 

More generally for M. 9 n one thinks of meromorphic differentials dQ A on T g of two 
types: (A) dVL a i is holomorphic outside P a with dVL a i = d[k % a + Oik' 1 )) near P a and 
(B) Q a (a 1) is a differential with simple poles at Pi and P a having residues ±1 
respectively so that 

dQ a , = dk a {k~ l + Oik' 1 )) = -dhik^ 1 + 0{k^)) (3.58) 

Here Q a ^(P) = J P dVL a ,i and for simplification of formulas later one complexifies 
the Whitham hierarchy leading to the moduli space M* gN = {T g , P a , fc" 1 , a^bi G 
Hi(T g ,Z)} where canonical homology basis. Then one normalizes the 

£Ia via § a . dfl A = and sets p(P) = fii i i(P) = f p dQiA- If one works on M. g N 
the normalization must be changed to £y § c dD, A = for c G Hi(T g ,7i). Now the 
multivalued function p(P) can be used as a coordinate everywhere on T g except 
at points U s where dp(U s ) = 0. One can accordingly change P to p in formulas 
such as ( |3.55| ) - ( |3 . 5 7] ) and we will do this without additional fuss. Thus a full 
system of local coordinates on M* N is given by {p a = p(P a ), v a ^ S) a = 1, ■ ■ ■ , N, s = 
— 1, 0, 1, • • • ; 7r s = p(Jl s ), s = 1, • • • , 2g, Uf = § b . dp, i = 1, • • • , g}. The associated 
compatible system of evolution equations will be 

d A k a (p,T) = {k a (p,T),n A (p,T)} (3.59) 

d A uf = d x uf (uf = <f dn A ); d A 7i s = d AP (u s ) = d x n A (u s ) (3.eo) 

J hi 



and we will say more about this below. 
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3.8 Zero curvature equations 

Now we have seen in Section 3.4 and in Section 3.6 how differentials of the form 
dp ~ Q\, dE ~ f2 2 , dfl ~ Q% arise in Whitham equations such as ( |3.54|) (or (|3.17|) ). 
Evidently further variables Tj could be appended with say Tj ~ Qi ~ k\ (and one 
expects fl at i will be associated to some T a ti as indicated in (|3.59| ) - ( p.60| ) - all of 



the relevant times are discussed below and the context of Hurwitz spaces is spelled 
out following f3T|). In any event one can expect to have a family of differentials 
Q A (A ~ (at,i) for example), such that Oa^b = 9 b Qa when the Qa are expressed in 
suitable variables. In connection with equations such as ( |3.59|) - (|3.60| ) (and (|3.17| ), 
(EJ751)) we make some observations on the possible forms of Whitham equations. In 
dealing with dispersionless KP (= dKP) for example there are two natural algebraic 
forms which emerge (cf. [|, [15|, |68], [74], f57|) which will be displayed below; further 
there are various geometrical developments in connection with Frobenius manifolds 



etc. a la [31], |32j . A (zero curvature) form which arises in this way is 



d A ^B - d B Vt A + {^a, ^b} = 0; {/, g} = f x g p - f P gx 



(3.61) 



where p ~ Q Ao {T Ao ~ X, A ~ (1, 1), p = Q^P) = / P dfii,i). We note that {f,g} 
here is {g, /} in |16|, |97[1 . Such equations ( |3.61|) can be regarded as compatibility 
equations for 

d A E = {E,n A ] (3.62) 

where E is an arbitrary function of (p, T) (easy exercise, using the Jacobi iden- 
tity). When d p E ^ one can write p = p(E,T) and d A f(p,T) = d A F(E,T) + 
(dF/dE)d A E (F(E,T) = f(p,T)). Then for such an E flgjnp becomes 



Indeed we note that (fl' A 



d A Q B (E,T) = d B Q A (E,T) 
(dQ A /dE), E' ~ (dE/dp)) 



(3.63) 



(3.6i) = o A n B + n' B E A - d B n A - n' A E B + (d x n A + n' A E x )n' B E' - (3.64) 
-n' A E'(d x n B + n' B E x ) = d A n B - d B n A + 

+n' B (E A + E'd x tt A ) - Q! A {E B + E'd x n B ) = 

But ( FTB2D implies for example that £4 = E x d p Q A {p, T)-E'd x Q A {p, T) = E X Q' A E'- 
E'(d x n A {E,T) + n' A E x ) = -E'd x tt A =>- E A + E'd x Q A (E, T) = 0. Similarly E B + 
E'd x Q B (E, T) = and (|3.64j) implies (|3.63| ). From this point, in coordinates E, T we 
can introduce a "potential" S(E,T) via (***) Q A (E, T) = d A S(E, T) (S ~ J2 Q A T A ) 
and write for Q = OS/ d B {8 ~ full exterior derivative) 



to 



SS{E, T) - QdE = d E SdE + £ d A SdT A - QdE = J2 ^ A dT A 



(3.65) 
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We note also that p = dSJ dX = Qa - It follows that 5uj = A &T& formally and 

from (|3.65|) one obtains 5uj = dE A 5Q ~ 5E A SQ. This indicates a very special role 
for E and Q = dS/ dE in the theory and this will be developed in a general way for 
dKP later following |L6j, |97| . In fact, by writing out 5u> (with = p), one will now 
have a "classical" string equation {Q, E} = 1 when Q, E are considered as functions 
of p, X (also then 8aQ = {Q, &a})- 



3.9 Times 

It is appropriate here to make a few comments about times Ta in a general sense (cf. 
here [||, 15, 17, 31]. Thus we will often have situations where say 



p 



(3.66) 



describe inverse functions. Often this arises in polynomial Landau- Ginzburg (LG) 
situations where e.g. p n+1 + a\p n ~ 1 + • • - + a n ~ A™ +1 and one thinks of Puiseaux series 
etc. but ( |3.66| ) can be more general. Then one notes that dp = d\(l — J2T mp m ^ 
so 



-m— 1^ 



Pr, 



Res p=00 X m dp 



.-Res x=00 X m ^-d\ 
m d\ 



(3.67) 



These coefficients p m will turn out to be universal time coordinates in a sense to be 
indicated below. • 



Let us consider now some examples from |68|, [74|] before bringing in the dKP back- 
ground. Thus one defines algebraic orbits of the Whitham hierarchy to be solutions 
Qa say obtained via a global solution of ( |3.62| ). For genus zero global means E is a 
meromorphic solution of (|3.62|) such that {E(p, T), k a (p, T)} = 0. 



3.10 Lax reduction, N = 1, g = 

Let Pi ~ oo with local parameter ki(p) = k(p) = p + J2T v s p~ s (recall p ~ P here). 
Suppose some power k n {p) ~ X n is a polynomial 

E = k n {p) =p n + u n . 2 p n - 2 + --- + u (3.68) 

We think here of k(p,T), E(p,T), etc. and it is convenient to take E = Ylxip — Pk) 
with distinct real roots pk leading to an n-sheeted Riemann surface T, E of genus 
associated with E (cf. Sections 4,5 and [IR pi). This means that the v s are 
(polynomial) functions of the Mj. The dispersionless Lax equations corresponding to 
this example will be, as in ( |3.62| ) 

diE(p,T) = {E(p,T),ni(j>,T)y, Qi(p,T) = [E$(p,T)) + (3.69) 
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(see also Section 4 for a more thorough development). We note here that from 
dip(E) = djc^i(E) as in ( ggg ) and (|3~69l) there results (using d x ^i(p) = d x ^i(E) + 
d E ni(E)E x and d^ip) = <M(£)£') 

d t E = Exd^i - E'(d x Qi{E) + d E n i (E)E x ) = -E'd iP (E) (3.70) 



Such an equation is also used in [|74|, ^5] but it's origin there is not clear (whereas 
here just its meaning is unclear - a shift X — > —X seems indicated). One defines now 
a generating function 

00 oo 

S = ]T TA = ]T Tik* + 0(/c _1 ); a = K (3.71) 

1 i 

Let q s satisfy (dE/dp)(q s ) = (i.e. (dE/dp) = np a ~ x + (n - 2)u n _ 2 p n ~ 3 + h 

U\ = Hi(p — q s ); q n = —J2i~Qs) an d consider q s , u as a new set of independent 
coordinates. Now let the Uj and Tj be related via (dS/dp)(q s ) = (this is a stipulation 
which determines a dependence of the Uj and Tj, leading to Uj = Uj(Tij). Then one 
shows that diS(E,T) = Qi(E,T) (the argument in |7|, [7S| is however somewhat 
curious and matters will be clarified in Section 4). We note also in passing from 
( p.69| ), namely diE = d x Ed p Qi — d p Ed x Qi, that at q s one has 

diE(q s ,T) = E x d p ni(q a ,T) (3.72) 

and this (Riemann invariant) form of (|3.69|) will be important later in topological 
field theory (cf. ||). Another point of view (related to Section 3.9) involves looking 
at any formal series Q(p) = Y^T^'jP' an d defining "times" via 

fi = -Mes p=00 k- i {p)Q{p)dE{p) (3.73) 

Thus for Q ~ (dS/dE) one has 

1 1 00 

Ti = -Res k=00 k~ l dS = -Res k=00 J2j T 3 kJ = T * ( 3 - 74 ) 
i i i 

Variations of this approach will appear in Section 4. It is clear however that one can 
use (|3.73| ) to define times Tj and as indicated in (|3.74p Q ~ dS/dE yields T = Tj. 
This leads to the more general picture below. • 

Take now the general situation with a "big phase space" as the moduli space 
Af g = {T g , dQ, dE, P a , n a , k' 1 , a^bi G Hi(T g )} where dE is a fixed normalized 
meromorphic differential having poles of orders n a + 1 at P a and dQ is a fixed nor- 
malized differential holomorphic outside of the punctures. Here E will have the form 

N n a 

E = p n + u n _ 2 p n - 2 + • • • + u + E ^s(P - Pa)' s (3.75) 

a=2 1 
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Generally in the Hurwitz space context E determines the ramification (cf. Section 
5.3). One also takes k™ a (p) = E(p) and dS = QdE. Then times are defined via 

T a>i = -Res Pa (k~\p)Q(p)dE(p)), (i > 0); T afi = Res Pa QdE (3.76) 

% 

and the formulas (i — 1, • ■ • , g) 



T hti = i dS- Tg j = -& dE- T E>i = i dQ (3.77) 

J cii J bi J bi 

Now one associates differentials to the time variables Ta- One should have as in ( |3.55| ) 
- ( 13.561) , Qa corresponding to Ta with § a . dVtA = 0, to which one adds holomorphic 



differentials dVLh,k normalized via 

d£l h ,k = $ik (3.78) 



to give the correct number of Ta (cf. |$2], |73|). The differentials dflE,i and dflQ^ are 
holomorphic on T g except for the dj cycles where they have jumps 

dVL% fi - dVL~ Ei = SijdE; dfij^ - dSl^ = 5ydQ (3.79) 

The normalization conditions are 

dn EJ = I dVL Qd = (3.80) 



Let then V g C M g be the subset of the big phase space where dE, dQ have no common 
zeros on T g . Then the system of times Ta above defines a system of coordinates on 
V g and the corresponding dependence of T g (T) and dE{T) on T represents a solution 
of the universal Whitham hierarchy on M* gN - This is equivalent to saying that the 
Whitham hierarchy can be considered as a way to define the special coordinate system 
on the moduli space of curves with punctures and jets of local coordinates in the 
neighborhoods of the punctures. The n a jets are the equivalence classes of coordinates 
where k' a (P) = k a (P) means k' a (P) = k a (P) + O^P)-™"- 1 ). 



3.11 The tau function 

We will make various remarks about the tau function here in the moduli space context 
and refer to Section 4 for other points of view. Generally the tau function of e.g. 
KP theory does not itself tend to a limit in dKP, but writing (in a WKB spirit) 
t = exp[(l/e 2 )F(T)} one can invoke limiting procedures and F(T) = logr dKP is the 
quantity of interest in the dispersionless theory (cf. || || |1^, [LJ| [18], |97| for discussion). 
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We follow here mainly |)S, |73|, [75[] (cf. also 32]). For genus zero and general E 
one can express solutions of (|3.62j ) implicitly via 



^(q s ,T) = 0; S(p,T) = J2(T A -T° A )Q A (p,T) 
dp Y 



(3.81] 



where (dE/dp)(q s ) = and for £> ~ (/3, z), dBS(p,T) = Qb(p,T). One defines then 

1 JV oo 



F = ~Y.( Res *{Y, f «AdS(p,T)} +T afi s a (T)) 
1 i i=i 



(3.82) 



where T Qj j ~ T a i — T° ^ and s a is defined via 



AT oo 



S(p,T) = £(£ + T afi logk a + s a ) + 0(k~ l ) 



(3.83) 



1 i=l 



(T 1;0 = — Y,2 T at o). One makes cuts now connecting Pi ~ oo with the P a and chooses 
a branch of S(p,T). Then 



s a = ^— (f log{p - p a )dS 
2m Jar, 



(3.84) 



where a a is a contour surrounding the corresponding cut (we will be somewhat cav- 
alier about interchanging p a and P a ). The function 5" has jumps on the cuts and is 
holomorphic outside of the cuts and punctures. One can also write fl3.82j ) in the form 

F = J dS AdS (3.85) 

which will be clarified as we go along (many details are missing in this first sketch a 
la Krichever). 

Now note from ( 3.82j ) for A = (a, % > 0) 

N oo 

2d A F = Res a (kidS) + £(£ Resp(T^dtt A ) + f^ A { P p)) (3.86) 

1 3=1 

where d A sp = VL A {pp). Now J2i Res a (Q A dQ,B) = implies that Resp(kj 3 d^l a)i ) = 
Res a (k z a dQ/3j) for j > 0. Further one obtains VL A {pp) — Resp(fi, A dlog(p — pp)) = 
Res a (k l a dVLf3fl). Putting these equations into (|3.86|) yields 



d aA F{T) = Res(k l a dS(p,T)) (i > 0) 



(3.87) 
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and similarly one obtains (4*) d a ^F(T) = s a . Hence the expansion of S(p,T) at P a 
has the form 



N oo 



S = E(E T ^ k i + T »,olog K + d a , F + E -d a>j F k; 



1 i=l 



which implies 



d% B F = ReScik^dSlB) {A =(a,i> 0)); 



(3i 



(3.89) 



d a ,odf3,oF = log(p a - p p ) 

Next, from ( ET57D and (4), Q becomes 2F = J2T A d A F from which 2F B 
2 TaFab + -Pb and hence we have 



2F = E ^£ = E E ^TbF^b 



(3.90) 



The quantities F^b will play the role of two point (correlation) functions in TFT 
and are in many respects the fundamental objects of the theory (see especially 

1, p|, H H, 0, H, 0) 

We look next at genus g and M g etc. as above. Then the moduli space Af g is foli- 
ated via leaves determined by the periods of dE, namely V k = § bk dE (corresponding 
to —Tq^ in (|3.77|) ). The discussion in ||| [F5[ is customarily mysterious here but 
we write down the following formulas involving F (clarification is needed). Thus one 
takes 

F = F Q + ^ E T E,kEdS - jf T Kk dS + T h>k T Eik E k (3.91) 

where F is given by (|3.82|) with T Q i replaced by T Q i . The first integral in ( |3.91| ) is 
taken over the left side of the a k cycle and E k = E(P k ) where P k is the intersection 
point of the a k and b k cycles; the last term with E k apparently makes F dependent 
only on the homology class of cycles. A picture can be contrived of course but more 
detail would be helpful here. For this F it is claimed that ( |3.87|) and (4) hold along 
with 

d h , k F = ^-[T Etk E k - <f dS}; (3.92) 

2m Jb k 

d E , k F = ^-.l EdS; d Q , k F = -^[<f QdS — 2T Ek T hk \ 
2m Ja k 4m Ja k 

Additional derivatives are given by 

d (h,k),A F = ^—[ E k^(E,k),A + Qk$(Q,k),A ~ j 1 ^ 



(3.93) 



24 



d(E,k),A F = ^~i k Edn ^ d (Q,k)A F = Q^a - d A {T Etk T Kk )] 



Further one obtains from (|3.93|) the b-period matrix of normalized holomorphic dif- 



ferentials on T g via 



= -f. ( 3 - 94 ) 

(cf. also ||). 



3.12 Landau-Ginzburg (LG) theory 

Let us go now to a Landau-Ginzburg (LG) situation with polynomial E(p) as in 
Thus work with N = 1, g = 0, and S+(p, 71, • • • , T m ) = £J»TA(p), ^ = 
£+ /n , E = k n = p n + ■ ■ ■ + uq, and think of = QdE = QE'dp or dS/d£ = Q (cf. 
G3.74Q , etc.). For m = n + 1 one will have a situation (cf. |75| ) d p S+ = (bip + bo)E' or 
9eS^ + i = b\p + b and one chooses &i = 1/n, bo = for the standard A n _i type LG 
theory (B n+ i = p/n in JFJ| - note also d p S+ = when £" = as desired in (|3.72|) etc. 



and see here also Sections 4 and 5). The LG theory here involves a superpotential 
W and a ring of primary chiral fields (pi isomorphic to R = C[p]/ {W'(p) = 0} where 
nW(p, T , ■ ■ ■ , T„_ 2 ) ~ E(p, T , 71/2, • • • , T n _ 2 /(n - 1), l/(n + 1)). Note the zeros of 
W correspond to zeros q s of E' so we are in the context of ( |3.72| ). With a change of 
variables Tj — > Tj_i/i then one asks for 

0, = -^-<9 p (nW0^ = -^c^+i = ±<W = ± J-0 i+1 £ (3.95) 
z + 1 i + l i + 1 

This is consistent with ( |3.72| ) again provided i?x = ±1 (which is OK at least in certain 
prototypical situations - cf. |], [17], ^| ) and indicates how the g s arise intrinsically 
since the ring R, as well as the time variables, are developed around the q s . Note 
however that E(q s ,T) ~ Riemann invariants and Ex is a natural creature in the 
Whitham equations (cf. ( |3.72| ) and ( p\4| )). In |3| one obtains a formula (cf. also 
0) 

d t F = Re Soo ( 2 n } " dp) (3.96) 
(z + l)(n + i + 



and from this it can be shown that 



-Re Soo (S- +1 dS+ +1 ) (3.97) 



where k n = nW and S n+1 = (k n+1 /(n + 1)) + Eo' 2 ( T j/U + After the change 

of variables 7} — > (Tj-i/j), S^ +1 agrees with the S^ +1 discussed at the beginning of 
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Remark 3.7. Now from <9e<S^~ + i(A;) = p/n one can write (|3.96| ) as 

9iF = Re SoD (— TT — — Trnd^gi) = (3 . 98) 
{i + l)(n + % + 1) dk n 

h.i+1 L.i+1 

= -Res 00 (—-d E S+ +l ) = -Res <x (r^-rdESn +1 ) 

l+l 2 + 1 

(note d(/g) = (df)g + /dp and Resd(fg) = 0). Let the right side of (|3.98|) be 
called Fx and compute diFi as follows. Using k(p) as local parameter one knows that 
diS+ +1 (k) = k^Ki + 1) = Qi/(i + 1) (since S+ +1 = m and % -> (T^/i, 

etc.). Hence diS n +i = k t+1 /(i + 1) is natural with diS~ +1 = k 1 ^ 1 /{i + 1). This gives 

m = -—^—Res^dS^ - k^dS- +1 ) (3.99) 

But = Res{k l+l dS n+ i) since dES n+ i = dk«S n+ i = (l/n)k~ n+1 dkS n+ i =>- k l+1 dES n+ i - 
(l/n)k l - n+2 d k S n+1 = d(J2a m k m ) and hence = Res^ 1 dS~ +1 + k^ l dS+ +l ). This 
gives then diFi = diF from (|3.98) and one will conclude that F — F\. 



4 CLASSICAL KP AND dKP 



We follow here |T3[ (cf. also |9"7| ) and begin with two pseudo differential 

operators (<9 = d/dx), 

00 oo 

L = a + ^n n+1 9- n ; W = l + Y,™nd- n , (4.1) 

1 i 

called the Lax operator and gauge operator respectively, where the generalized Leib- 
nitz rule with d~ l d = dd~ l = 1 applies 

<y.f xf'V/^' 1 (4.2) 



3=0 



for any j £ Z, and L = WdW 1 . The KP hierarchy then is determined by the Lax 
equations (d n = d/dt n ), 

d n L = [B n , L] = B n L - LB n , (4.3) 

where B n = L" is the differential part of L n = L% + L n _ = t%& + ElL C^- One 
can also express this via the Sato equation, 

dnWW- 1 = -L n _ (4.4) 
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which is particularly well adapted to the dKP theory. Now define the wave function 
via 

CO CO 

4, = Wet = w(t,\)et; e = E^ A "; w(t,X) = l + Y / Mt)^~ n , (4-5) 

1 1 

where t\ = x. There is also an adjoint wave function if)* = W* _1 exp(— £) = 
w*(t, A) exp(— £), w*(t, A) = 1 + J2t° w*(£)A~\ an d one has equations 

Lip = XiP; d n ifj = B n iP; = Xi/j*; d n iP* = -B* n iP* . (4.6) 



Note that the KP hierarchy (|4.3|) is then given by the compatibility conditions among 
these equations, treating A as a constant. Next one has the fundamental tau function 
r(t) and vertex operators X, X* satisfying 

, (t , A) = *»m = = fn^pi; (4.7, 

T{t) T{t) T{t) 

ib*(t X) = X *W T ® = e ' €g +(A)r(t) = e'Mt+lA- 1 ]) 
1 ' ) r(t) r(t) r(t) 

where G±(X) = exp(±^(9, A" 1 )) with 8 = (d u (l/2)«9 2 , (l/3)<9 3 , • • •) and t ± [A" 1 ] = 
(ti ± A _1 ,t 2 ± (1/2)A -2 , • • •). One writes also 

COO \ CO 

XX A ") =EXj(M* (4.8) 

where the Xj are the elementary Schur polynomials, which arise in many important 
formulas (cf. below). 

We mention now the famous bilinear identity which generates the entire KP hier- 
archy. This has the form 

r ip(t,X)ip*(t',X)dX = (4.9) 



where ^(^dX is the residue integral about 00, which we also denote Res\[(-)dX]. 
Using (|4.7j ) this can also be written in terms of tau functions as 

r(t - [A- 1 ])^' + [X-^e^-^'^dX = (4.10) 

This leads to the characterization of the tau function in bilinear form expressed via 
(t _ t - y, t' -> t + y) 

(jr X n(-2y)Xn + i(d)e^yA r ■ r = (4.11) 
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where dfa ■ b = (d m /dsf)a(tj + Sj -)6(tj- - s^Uo and 9 = (ft, (l/2)<9 2 , (1/3)3,, • • ■)• 
In particular, we have from the coefficients of y n in (|4.11 ), 

dAr ■ r = 2 X n + i0)r ■ r (4.12) 

which are called the Hirota bilinear equations. One has also the Fay identity via (cf. 



T5| - c.p. means cyclic permutations) 

E( s o - si)(s 2 - s 3 )r(t + [s ] + [si])r(t + [s 2 ] + [s 3 ]) = 



(4.13) 



c.p. 



which can be derived from the bilinear identity ( 4.10|) . Differentiating this in sq, then 
setting s = s 3 = 0, then dividing by S]S 2 , and finally shifting t — > t — [s 2 ], leads to 
the differential Fay identity, 

T(t)dr(t + - [s 2 ]) - r(t + - [s 2 ])dr(t) 

= (^r 1 - s, 1 ) [r(t + [ 8l ] - [s 2 ])r(t) - r(t + [ Sl ])r(t - [s 2 ])} (4.14) 

The Hirota equations ( |4.12| ) can be also derived from ( 4.14Q by taking the limit 
s% — > s 2 . The identity ( 4.14 ) will play an important role later. 

Now for the dispersionless theory (dKP) one can think of fast and slow variables, 
etc., or averaging procedures, but simply one takes t n — > et n = T n (ti = x — >• ex = X) 
in the KP equation u t = {1/A)u xxx + 3uu x + (3/4)9 _1 m to , (y — t 2 , t — t 3 ), with 
d n -> e<9/<9T n and w(t n ) -> f/(T n ) to obtain <9 T ?7 = 3f/f/ x + (3/A)d- 1 U YY when 
e — > (d = d/dX now). Thus the dispersion term u xxx is removed. In terms of 
hierarchies we write 



L e = ed + Y,Un + i(T/e)(edY 



(4.15) 



and think of u n (T/e) = U n (T) + 0(e), etc. One takes then a WKB form for the wave 
function with the action S 

■1 



if) = exp 



-S(T,X) 



(4.16) 



Replacing now d n by ed n , where d n = d/dT n now, we define P = dS = Sx- Then 
€ l d l if) — > P l ip as e — > and the equation = becomes 



(4.17) 



where the second equation is simply the inversion of the first. We also note from 
d n ip = B n ip = J2q b nm (^d) m ip that one obtains d n S = B n (P) = A™ where the subscript 
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(+) refers now to powers of P (note ed n ip/ip — > d n S). Thus B„ = L" — > B n (P) 
A+ = So b nm P m and the KP hierarchy goes to 



(4.18) 



which is the dKP hierarchy (note = B n => d„P = dE n ). The action S in ( [3.13D 
can be computed from ( |3.5| ) in the limit e — > as 



(4.19) 



where the function F = F(T) (free energy) is defined by 

1 



t = exp 



;F(T) 



(4.20) 



The formula ( |4.19| ) then solves the dKP hierarchy fl4.18|) , i.e. P = B\ = dS and 



B n = d n S = X n ~J2 



711 



(4.21) 



where F nm = d n d m F which play an important role in the theory of dKP. 



Now following [97| one writes the differential Fay identity (^4.14f) with ed n replacing 
d n , looks at logarithms, and passes e — > (using fl4.2U| )). Then only the second order 
derivatives survive, and one gets the dispersionless differential Fay identity 



E ^~ mx ~ 



m,n=l 



mn 



l0 T"? 



n 



(4.22) 



Although ( [4.22|) only uses a subset of the Pliicker relations defining the KP hierarchy 
it was shown in |)7| that this subset is sufficient to determine KP; hence Q4.22D 
characterizes the function F for dKP. Following [15, 18 1, we now derive a dispersionless 
limit of the Hirota bilinear equations (|4.12| ), which we call the dispersionless Hirota 
equations. We first note from ( |4.19| ) and ( ^4.17| ) that F ln = nP n+ i so 



-n F lr 



n 



E^+iA- n = A-P(A) 



Consequently the right side of ( |4.22| ) becomes logf ^ _^ ] and for fi 
P = d\P we have 



(4.23) 
A with 



logP(A)= E A ~ 



m,n=l 



■m-n ^ "in 

mn 



E E 

j=l \n+m=j 



\~ 3 



mn 



(4.24) 
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Then using the elementary Schur polynomial defined in Q4.8| ) and ( |4.17| ), we obtain 



P(X) = £ Xj (Z2, • • • , Z s )\-i = 1 + E Fi^- 1 ; 
o 1 



E 



mn 



(Z 1 = 0) 



m+n=i 

Thus we obtain the dispersionless Hirota equations, 



F\j — Xj+i{Z\ — 0, Z 2 , 



(4.25) 



(4.26) 



These can be also derived directly from (|4.12|) with (|4.20 ) in the limit e — > or 
by expanding (|4.24j) in powers of X~ n as in JI5], [18|]). The equations (|4.26|) then 
characterize dKP. 



It is also interesting to note that the dispersionless Hirota equations fl4.26|) can be 
regarded as algebraic equations for "symbols" F mn , which are defined via ( |4.21| ), i.e. 



B„ 



K 



a"-E 



nm x -m 



m 



(4.27) 



and in fact 



Resp[X 



m d\ n + ] 



(4.28) 



Thus for A, P given algebraically as in (|4.17|) , with no a priori connection to dKP, 
and for B n defined as in ( 4.27| ) via a formal collection of symbols with two indices F mn , 
it follows that the dispersionless Hirota equations (|4.26|) are nothing but polynomial 
identities among F mn . In particular one has from |JT5| 



THEOREM 4.1. 

hierarchy. 



.28|) with ( |4.26| ) completely characterizes and solves the dKP 



since 



Now one very natural way of developing dKP begins with ( |4.17| ) and (^4.18|) 
eventually the Pj+\ can serve as universal coordinates (cf. here for a discussion 
of this in connection with topological field theory = TFT). This point of view is also 
natural in terms of developing a Hamilton- Jacobi theory involving ideas from the 
hodograph — Riemann invariant approach (cf. [IB, |6, |62], B31 and in connecting 



NKdV ideas to TFT, strings, and quantum gravity. It is natural here to work with 
Q n := {l/n)B n and note that d n S = B n corresponds to d n P = dB n = ndQ n . In this 
connection one often uses different time variables, say T' n = nT n , so that d' n P = dQ n , 
and G mn = F mn /mn is used in place of F mn . Here however we will retain the T n 
notation with d n S = nQ n and d n P = ndQ n since one will be connecting a number 
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of formulas to standard KP notation. Now given ( |4.17| ) and ( |4.18| ) the equation 
d n P = ndQ n corresponds to Benney's moment equations and is equivalent to a system 
of Hamiltonian equations defining the dKP hierarchy (cf. fTE| , |62|); the Hamilton- 
Jacobi equations are d n S = nQ n with Hamiltonians nQ n (X, P = dS)). There is now 
an important formula involving the functions Q n from ||62|| , namely the generating 
function of dpQ n (X) is given by 



In particular one notes 



POO - P(A) 



J2d P Q n (X)fi- n 



(4.29) 



u' 1 

P{») - P(A) 



djM = d P Q n+ i(\) 



(4.30) 



which gives a key formula in the Hamilton- Jacobi method for the dKP [ W2\ . Also note 
here that the function P(A) alone provides all the information necessary for the dKP 
theory. It is proved in [jTHj that 



THEOREM 4.2. The kernel formula (|4.29| ) is equivalent to the dispersionless 
differential Fay identity ( |4.22|) . 



The proof uses 



dpQn — Xn-l(Ql, ' ' " ) Qn- 



where Xn(Qi, ■ • • , Qn) can be expressed as a polynomial in Qi 
cients given by polynomials in the Pj+i- Indeed 



(4.31) 

P with the coeffi- 



Xr. 



det 



P -1 

P 2 P 

P 3 P 2 

Pn Pn-1 







p 




















Pa Ps P 2 P 



d P Q 



71+1 



(4.32) 



and this leads to the observation that the F mn can be expressed as polynomials 
in Pj + i = Fij/j- Thus the dispersionless Hirota equations can be solved totally 
algebraically via F mn = $ mn (P 2 , P 3 , • • • , P m+n ) where $ mn is a polynomial in the 
Pj + i so the Fi n = nP n +i are generating elements for the F mn , and serve as universal 
coordinates. Indeed formulas such as (f4.32|) and ( |4.31|) indicate that in fact dKP 
theory can be characterized using only elementary Schur polynomials since these 
provide all the information necessary for the kernel (f4.29|) or equivalently for the 



dispersionless differential Fay identity. This amounts also to observing that in the 
passage from KP to dKP only certain Schur polynomials survive the limiting process 
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e — ► 0. Such terms involve second derivatives of F and these may be characterized 
in terms of Young diagrams with only vertical or horizontal boxes. This is also 
related to the explicit form of the hodograph transformation where one needs only 
dpQn = Xn-i(Qi, 4 ■ ■ , Qn-i) and the Pj+\ in the expansion of P (cf. JT5J). Given KP 
and dKP theory we can now discuss nKdV or dnKdV easily although many special 
aspects of nKdV for example are not visible in KP. In particular for the one will 
have F nj = F jn = for dnKdV. We note also (cf. H) that from ( pL3l ) one has 

1 oo oo oo 

p , v pm = E d P Qn^ n = E Xn(Q)^ n = exp(E Qm^ m ) (4.33) 

- ^\ A ) 1 1 

5 HURWITZ SPACES 

The preceeding development leads one to work in the framework of Hurwitz spaces and 
Frobenius manifolds developed by Dubrovin (cf. |TLf). In fact we could have started 
here by hindsight but felt it desirable to pursue the more tortuous path traced in 
order to illustrate various points of view, derivations, and techniques (in particular 
the nature of averaging) . We will extract here freely from |31[ . 



Hurwitz spaces are defined to be moduli spaces of pairs (E g ,A) where E g is a 
smooth algebraic curve of genus g and n + 1 sheets; A is a meromorphic function on 
E g of degree n+1 (which can be used to realize E g as an n-sheeted covering over CP 1 ). 
One will assume here that the ramification over oo is fixed as indicated below. Thus 
let M = M 9;n0j ... „ m be a moduli space of dimension n = 2g + n + - ■ ■ + n m + 2m of sets 
(E 9 ; oo , ■ • • , oo m ; A) where E g is a Riemann surface with marked points oo , ■ ■ • , oo m 
and A is a meromorphic function A : S 9 —>■ CP 1 with A _1 (oo) = oo U ■ ■ • U oo m and 
having degree rii + 1 near oo$. The critical values of A, defined via 

u j = \(Pj); dX\p. = 0(j = l r --,n) (5.1) 

are to be local coordinates in open domains M where ^ u % for i ^ j (here M is to 
denote a covering of our eventual Hurwitz space in which A is a variable). These are 
the ramification points of the surface A : E — > CP 1 and the Pj are branch points of 
S g . One asks also that the one dimensional affine group acts on M via 

(E g ] oo , ■ ■ ■ , oo m ; A) -> (E g ; oo , • • ■ , oo m ; aX + 6); u l i-> au l + b (5.2) 

We note that X(Pj) = is a point in a moduli space while Pj refers to a particular 
Riemann surface. There is a strong interaction between E g and the function A as 
indicated in the examples to follow. 
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5.1 Examples 

(A) For g = 0, m = 0, uq = n the Hurwitz space consists of all polynomials of 
the form A(p) = p n+1 + a n p n ~ x + • • • + oi (flj 6 C) (here p denotes a point on E g 
and strictly one should use P - cf. below). The affine transformations A i— > aX + b 
act via p i— > a 1 ^ n+1 - ) p, cij i— > aja( n-i+1 )/ (n+2 ) (i > 1), ax i— > aai + &. This corre- 
sponds to A n as in Section 5.3 and the extension to M gn is displayed in Section 
7.3. (B) For g = 0, m = n, n Q = ■ ■ ■ = n m = the Hurwitz space is all ra- 
tional functions of the form X(p) = p + Yli[li/(p ~ Pi)] with affine group action 
p i— > ap + b, pi f— > + (b/a), i— > ao;. (C) For p > 0, m = 0, n = 1 the Hur- 
witz space consists of all hyperelliptic curves fi 2 = ni 9+1 (A — u^). The critical values 
it 1 , - • • , u 2s+1 of the projection (A,/i) ^ A are the local coordinates on the moduli 
space (note /z 2 = Poly(X,u j ) = P(X,u j ) =>- 2ji = d\P(dX/ 'dp, =>■ dX/dpi = for fi = 
or A = it-') . Other examples are indicated below. • 

Evidently example (C) is defined somewhat differently and we will try to clarify 
this later. One notes here from |59| for example that if S g is the Riemann surface of 
an irreducible algebraic equation P(z, w) = of degree n'mw and if the branch points 
have orders rii then the genus is p = 1 — n + (1/2) Y7\ ^« (Riemann-Hurwitz formula). 
Thus e.g. for A(p) in the first example (A) (with p ~ w) if one assumes A'(p) = 
Eli (p — Ps) an d distinct p s then n s — 1 and there is a branch point of order n at oo. 
Hence p = 1 — (n + 1) + (1/2) 1 + (l/2)n = as required here. This shows also how 
an n-sheeted surface can have genus 0. For hyperelliptic situations we recall the fol- 
lowing facts (cf. ||102|| ). Consider R = fi 2 = ]lo m (^ — ^fc) as i n the beginning of Section 
5.1 with branch points Xk and oo. One makes cuts (Ai, A2), (A3, A4), • ■ ■ , (\2m-1, X2 m ) 
as spectral gaps with Oj cycles corresponding to (A2i-i, A2i) around these cuts verti- 
cally and bj cycles horizontal around (Aq, Ai), • • ■ , (A2 TO > 00 ) f° r example (see 0, |102|| ). 



This produces a surface with p = m visible holes. Note the Riemann-Hurwitz formula 
gives also (w ~ H) g = 1 - 2 + (1/2) Eo m 1 + (V 2 ) = m - Th e hyperelliptic situation 
will be developed further below. 

Now one constructs a covering M of M = M 9;n0i ... jnm in order to describe multi- 
valued quadratic differentials for which the one forms Qq = J2i du % Res p^Q / dX) will 
define metrics on the Hurwitz space via Q(a ■ b) =< a,b >n- Thus take M to con- 
sist of sets (E s ; ooq, ■ ■ ■ , oo m ; A; ko, ■ ■ - , fc m ; a\, ■ ■ ■ , a g , b\, ■ ■ ■ , b g ) G M with the same 
(E g , oo , ■ • • , oo m , A) as before and with a marked symplectic basis dj, fej G Z) 
and marked branches of roots near oOj of orders rij + 1 (i = 0,---,m), namely 
fc™ i+1 (P) = A(P) near ooj. The admissable quadratic differentials will be constructed 
as Q = <j) 2 for certain primary differentials cf) on S g or on a covering. The primary 
differentials have the following forms. (D) Normalized Abelian differentials of the 
second kind on E g with poles only at ooq, • • • , oo m of orders less than the correspond- 
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ing orders of the differential dX. Explicitly = <p t i-,a (i — 0, - • • , m, a — 1, • • • , nj) is 
the normalized Abelian differential of second kind with a pole at ooj of the form 



dk^ + regular terms near oof, (f <p t i; a = (5-3) 

a Ja,j 



fE) Next one considers 



^2di<f) v i (i = 1, • • ■ ,m) (5.4) 



with the 5i independent of the point in M. Here 0^ is one of the normalized Abelian 
differentials of second kind on S 9 with a pole only at ooj with principal part of the 
form 

(j) v i = —d\ + regular terms near oc>j; ® cf) v i = (5-5) 

Jaj 

(F) Next take 

m „ 

= J2ai(j) w i] f = (5.6) 

1 Ja i 

with ai, • ■ • , a m independent of the point on M and <f) w i is the normalized Abelian 
differential of the third kind with simple poles at oo and ooj with residues —1 and 
+ 1 respectively. (G) Next look at = J2i/3i4>r i with /3j independent of the point in 
M and r i is the normalized multivalued differential with increments along the cycles 



(P + bj) - <j) ri {P) = -6ijd\; A r , = (5.7) 



bj of the form 



and without other singularities. (H) Finally = J]? li ( t ) s i with 7^ independent of 
the point on M and <p s i denotes holomorphic differentials normalized by § a . <ft s i = 5^. 
Then one will pick a primary differential Q = 2 to develop a Frobenius structure 
on M. In constructing the superpotential one introduces a multivalued function p on 
S 3 via p(P) = v.p. J^ o where divergent parts have been subtracted in the principal 
value integral. Then = dp and A = X(p) can be used. 

5.2 WDVV and Frobenius manifolds 

We prepare the way by extracting some material on the WDVV equations and Frobe- 
nius manifolds from [[H], 32|. We will sketch (without proofs) only a brief selection 



of this extensively developed material. The main creature is a quasihomogeneous 
function F(t), t ~ (ti, ...,t n ) such that the third derivatives d a d^d-yF = c a/ 3 7 (t) cor- 
respond to correlation functions < a 0/30 7 >o and in the LG models exp(F) is the 
dispersionless tau function (we use t^ and T k interchangably here - T k is correct in our 
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notation but tk is used in |3T], 0). These functions c a p 7 satisfy a system of partial 
differential equations (PDE) called the WDVV (Witten- Dijkgraaf-Verlinde-Verlinde) 
equations based on algebraic properties of TFT. Thus (cf. |H, one takes a 2-D 



TFT with say n primary fields (pi and < , >o ~ genus zero correlation function. The 
two point functions < 4> a 4>p >o= r\ a p = j]p a determine a nondegenerate scalar prod- 
uct (see below) on the space of primaries and the c a/ 3 7 determine the structure of the 
primary chiral algebra. More axiomatically one looks for F(t) with c a ^ = d a dpd 7 F 
such that ri a p = c\ a p is a constant nondegenerate matrix with rf^ = {7] a p)~ l . Then 
one specifies clg = rp^c^p and such c a p^{t) must define in the n-dimensional space 
with basis (1 < i < n) an associative algebra A t with 

e a -e p = cl p (t)e 7 ; cP la {t) = 8% (5.8) 



Here F should be quasihomogeneous in the sense F(c dl ti, ...,c dn t n ) = c dp F(ti, ...,t n ) 
(we also use tj and f interchangably) and this can be expressed infinitesimally via 
E = E a {t)d a and (£ ~ Lie derivative) 

C E F = E a (t)d a F(t) = d F ■ F(t); E a = d a t a (5.9) 

Note for e = d\ (~ unity vector field) = —die. One can modify quasihomogeneity 
by requiring 

C E F = dpF + A^rt 13 + B a t a + C (5.10) 

since the third derivatives are not affected. Further if dp ^ 0, dp — d a ^ 0, dp — 
d a — dp ^ for any a, (3 then the extra terms in ( |5.10| ) could be killed by adding a 



quadratic function to F. The degrees dx, ...,d n , dp are well defined up to a nonzero 
factor and we consider here only the case di / with normalization so that d\ — 1. 
Often one writes d a — 1 — q a (q± = 0), dp = 3 — d, so q n = d and q a + g„_ a +i = d. 
Then associativity implies the WDVV equations 

dJdpdxFrj^d^dsd^F = d^dpdxFrj^d^sd^F (5.11) 

for 1 < a, f3, 7, 5 < n. A solution F of ( |5.11| ) is called a primary free energy. 



Now one says A is a (commutative) Frobenius algebra if (1) A is a commutative 
associative algebra with unity e (2) (a, b) \— >< a, b >: A x A — > C is C bilinear 
symmetric nondegenerate with < ab,c >=< a, be >. Then < a,b >=< e,ab > 
(e ~ ei) will be a suitable scalar product and for e« a basis in A one can specify 
< ei,ej >= rjij, eiej = c^eu to obtain a structure of Frobenius algebra on A (cijk = 
rji S c S j k = Cjik = Cikj). A theory of Frobenius manifolds can now be developed as follows. 
M is a Frobenius manifold if a structure of Frobenius algebra is specified on every 
tangent space T t M depending smoothly on t such that (1) < , > is a flat metric (2) 
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Ve = for the Levi-Civita connection V of < , > (3) For c(u, v, w) =< u- v,w > one 
requires (V z c)(u,v,w) to be symmetric in u,v,w,z and (4) V(V-E) = for a vector 
field E = Euler vector field (cf. ( 5.10Q ) such that the corresponding one parameter 



diffeomorphism group acts via conformal transformations on the metric < , > and by 
rescalings on the F-algebra T t M. The infinitesimal form of (4) is 

V 7 {V p E a ) = 0; C E cZp = c%- C E e = -e; C E r] a p = D Va p {D = 2 - d) (5.12) 

One proves that any solution of WDVV with d\ ^ defined in a domain of t deter- 
mines a structure of F-manifold via 

d 

d a -dp = clpdy, < d a , dp >= r] af3 ; d a = — ; e = d 1 (5.13) 

with E of ( [5.101 ). An example of particular interest arises from M = {X(p) = p n+1 + 
o>nP n ~ x + ■ ■ ■ + cli} with T\M ~ polynomials of degree < n. Then the Frobenius 
algebra A x on T X M is C[p}/(\'(p)) with < f,g > x = Res p=00 [f(p)g(p)/X'(p)}. Here 
e = d/dcti and E = (1/n + 1) J2( n — i + l)a>idi. 

5.3 Aji—\ 

We will write out first the situation for the topological minimal model following 



1 32 1 since it is more detailed than 0] and will give a more complete picture. There 



is a lot of duplication with previous formulas but we will stay here with the notation 



of [Q. The dispersionless theory will also be presented again in this notation. Thus 
one considers the coupling space M of all polynomials of degree n of the form M = 
{X(p) = p n +a n -2p n ~ 2 + - ■ -+ao} and the associated F-algebra is A\ = C\p]/(X'(p) = 0) 
with scalar product < f(p),g(p) >= —(^/n)Res p=00 [f(p)g(p)/X'(p)] (p ~ P). There 
is an affine structure on M defined as follows. Let <p a {p, X), a = 1, ...,n — 1, be an 
orthogonal basis of A\ satisfying 

< 4> a , 4>/3 >= Vaf3 = 5 a+f 3, n ; deg 4> a = a - 1 (5.14) 
Generally fields can be defined via 

Ti — 

(j> a (p,X) = -d p Xl (5.15) 
a 

The dependence of A(p) (i.e. of its coefficients) on the flat coordinates t a (better T a ) 
is determined via 

d a X(p) = -<p a (p, A); a= l,---,n- 1 (5.16) 
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We recall the Gelfand-Dickey hierarchies now via d q L = [L, L+] for L = d n + 
a n _ 2 <9 n ~ 2 + - • -+ao- The dispersionless hierarchy is formed as before via X = ex, T q = 
et q , e — > and one obtains in the spirit of |39| (d q ~ d/dT q now and d ~ d/dX) 

d a dp\x= c = d x d$ a \ x=c ; $ Q = A| (5.17) 

This can be rewritten via Wronskians as 

d a X(p) = d p \(d x $ a )\ p=k - d p $ a (d x \)\ p=k (5.18) 

(these are all standard calculations in the dispersionless theory (cf. [0). Now one 
has diagonal coordinates (Riemann invariants) for our dispersionless hierarchy as the 
critical values U\, ■ ■ ■ , u n -\ of X(p), i.e. 

Ui = \{ Pi ); A'(pi) = (5.19) 



and the characteristic speeds are (cf. also f3"2fl , equation (5.58), where d p & q ~ 
— (Q q /fli), Q\ ~ dp in the small phase space) 

v q ,i(u) = -dp® q \ p=Pi ; d q u l = v q $ x u l {% = 1, • • • , N = n - 1) (5.20) 

We assume here that A'(p) = has simple zeros so the algebra A\ will be decompos- 
able by definition and the corresponding metric ds 2 in these coordinates is 

(5 ' 21) 

Now one knows that the functions on M 

h a>q = --^-—Res p=00 X^ /n)+q dp; a = 1, - • • , n - 1, q > (5.22) 

where (a) q = a(a + 1) ■ ■ ■ (a + q — 1), are basic conservation densities and 

N Q 

dh a , q = h a , q -i (q > 1); dh afi = c (d = ^<9j, d { = — ) (5.23) 

(recall s n ~ Res L n ~ nH\__ 1 and Res L n — > Resp\ n so ( |5.22|) is natural enough). 
This follows from noting that a translation A i— > A + e, p i— >■ p, <2j i— > aj (i ^ 0), a i— > 
a + e is equivalent to translation along + • • • + u x so that 

<9/i aj g = D e h a>q (ui + e)| e=0 = (5.24) 
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—^—Res p D e [(\ + e)^ n ^dp(X + e)] e=0 = Res^^+^dp 

(a/n) q " (a/n) q -i 

This gives ( |5.23j ) for q > 1 and dh a> o = Res p X^ a ^ n ^~ l dp = —5 a ^ n -i. In particular for 
the flat coordinates we have 

T a = —nReSp dp; a = 1, ■ ■ ■ , N = n — 1 (5.25) 

n — a 

(cf. Remark 3.5 and Example 3.6). One has in fact a generating function for the h a q 
in the form 

xJt, z) = Res p iFAl; 1 + -; zX)dp; a = 1, • • • , N (5.26) 

a n 

where \F\{a\ c; z) = X)o° (°) m m! ^ s ^ ne Kummer hypergeometric function. 

Now returning to |yj let be one of our primary differentials and let be 
the open domain in M specified by the condition 4>(Pi) ^ (i = l,---,n). The 
main result says that for any primary differential <fi of the types (D) to (H) the 
multiplication <9j ■ dj = 5ijdi for <9j = d/du 1 , the unity and Euler vector field e = J2± 
and E = J2iU l di, and the one form Vt^ determine on a structure of Frobenius 
manifold. The flat coordinates T A , A = 1, • • • , n, consist of the five parts 

T A = ^ T i;a^ o < z < m, 1 < a < n; q j , 1 < j < m; r j , s j , 1 < j < g) (5.27) 
where 

T' a = Res^k^pdX, i = 0, • • • , m, a = 1, • • • , n { \ (5.28) 
p l = v.p. / dp, q l = -Res^.Xdp, i = 1, • • • , m; 

*/ oon 



r l = / dp, s l = -— / Adp, i = 1, • ■ • , g 
The metric will have the form 

7] T i ;aTJ ,f3 = — —— Sijd a+ /3 t n i+ i; (5.29) 

TLi -\- 1 

1 x . . _ J_r 

The function A = X(p) will be the superpotential of this Frobenius manifold in the 
sense that u l {T) = X(q l (T),T) where (dX/dp)(q l (T)) — (i — 1, ■■■,n), plus other 
stipulations which will be indicated as needed. Finally for any other primary differ- 
ential d> the one form is an admissable one form on the Frobenius manifold in a 



38 



sense to be indicated. 



We check some of this following |pl| . First one declares that the u l (1 < i < n) 
are the canonical coordinates for a multiplication di ■ dj = 5ijdi (di = d/du % ) and 
defines the metric corresponding to via < d',d" >0= ^(<9' ■ d") for any two 
tangent vector fields d', d" on M . This metric will also be diagonal in these canonical 
coordinates and one writes 

ds\ = Y^riii(du 1 ) 2 ; rj ti = Res Pt — (5.30) 

One can show now that, for any primary differential, this is a flat Darboux-Egoroff 
metric in an open domain leading to a Frobenius structure on M^. This holds for 
any <p and UM^ provides a Frobenius structure on M. The proof is nontrivial and we 
refer to W%- Next one shows 



d T i;c\(p)dp = — T * ;a ; (5.31) 

d v i\(p)dp = — (fi v i; d w z\(p)dp = —0^; 
d r i\{p)dp = —<fi r i; d s i\{p)dp = —<fi s i 
This is achieved via the thermodynamic identity 

d a (Xdp)\ p=c = -d a (pd\)\ x =c> (5.32) 



For (|5.32|) note that A = X(p(\,T),T) implies = d a X\ P = c + (dp / d\)d a p\\ =c i which 



says d a X\ p=c dp = —d a p\\ =c >d\ which is what (|5.32|) means. Now for example deriva- 
tives (dp(\) / dT A )\\ =c i are holomorphic on S 9 /oo away from the Pj where they have 
simple poles, and d\ vanishes precisely at such Pj so ( |5.32| ) is holomorphic on £ 9 /oo. 
A calculation of singularities gives ( |5.31|) . Finally let T A be one of the coordinates 
( |5.28| ). Then to prove ( |5.29| ) look at (f>A = —d T A\dp. Some calculation (not shown 
here) shows that 

< d T A,d T B > = Res dX=0 ~~7T~ = d e < (j)A<pB > (5.33) 

|A|<oo dX 

where < > is a complicated explicit bilinear pairing of suitable differentials, fur- 
ther nontrivial calculation yields ( [5.29|) . The structure constants for the Frobenius 
structure will now be 

" <j) T A<j) TB (f) T c 

c ABC = YRes Pi ^ (5.34) 
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and the Frobenius structure can be extended to the entire moduli space M via the 
condition that (</> t a</> t b — c ( ^ B <f) T cdp)/d\ be holomorphic for \w\ < oo. It turns out 
also that 

F — — - < pdX pdX >; d T Ad TB F = - < <p A <pB > (5.35) 

In particular 

d s ad s i3F = -T a p = f <p s c (5.36) 

This means that the Jacobians J(£ 9 ) = C 9 /m + rn, m, n e Z 9 are Lagrangian mani- 
folds for the symplectic structure J]? ds a Adz a where natural coordinates 
on J(Sp) coming from the linear coordinates in C 9 . Thus the Jacobians are complex 
Liouville tori and the coordinates z a , s a are complex action-angle variables on the 
tori. 



6 CAUCHY TYPE KERNELS 

6.1 Review of KdV 

For convenience we begin with a hyperelliptic Riemann surface (corresponding to the 
KdV situation) in the form R 2 (X) = Flo 9 (A - A fe ) as in Remark 3.3. Then there are 
normalized holomorphic differentials 

Ul ~ — rT\) — ' t Ui = 



(again u>i is used instead of duJi). We refer here to || |[ |7], [|, |17[ |39], |51j, [|^, ^3 



for hyperelliptic surfaces. There is often some difference in notation or normalization 
in the literature and we will try to specify explicitly our selection here (to follow 
102|| ). Thus with R 2 as above one has branch points A , • • • , A 29 , oo (whereas for 
R 2 = I\o g l (X— Afc) there are branch points A , • • • , A 29 _i). Again ~ (A 2i _i, A 2 «), i = 
1, • • • , g, and we will typically take Aj real and finite with Ao < Ai < • • • < A 2r The 
period matrix is Bij = uJi. We note that in |51], [|5[] for example one uses cuts 
(— oo, Ao), • ■ ■ , (A 29 _i, A 29 ) instead of (Ao, Ai), • ■ ■ , (A 29 , oo) as is indicated here; such 
choices are obviously equivalent. The 6« cycles can be drawn e.g. from a common 
vertex P passing through (A 2j _i, A 2j ). We recall that hyperelliptic Riemann surfaces 
have a number of very special properties (cf. || f4(| f|8], |53|, |84l ) so some statements 



to follow have correspondingly limited applicability. We keep our standard local 
coordinate k(P) or k~ l (P) at the single point P\ ~ oo. We use differentials Qk ~ ^i,fe 
as in Remark 3.4 with dVlk ~ d(k t + 0(k~ 1 )) near oo, § dVlk = 0, andp(P) = J p d£l\. 
Set then Uj = § h dflj as at the beginning of Section 3.1 (mod 2m - U\ ~ Uk, ~ 
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V k , Ujf ~ Wfe). Recall also U t A = § b% dfl A and A ~ (1, /c) here with dp = Q 1 , dE = 
Q2, dQ = Q3. 



Now in the hyperelliptic case the u>k have the form ( |6.1|) as indicated earlier. We 
observe that the Riemann surface is generated e.g. via the spectral bands of Bloch 
eigenfunctions as in Remark 3.3 and has nothing to do a priori with a putative LG 
potential as in ( |3.68| ) (cf. however (C) in Section 5.1). For g 7^ one expects tO{ as in 
( |6.1|) with atij determined by the normalization conditions § uji = 5ik (g equations 
in g unknowns for each i). Note also as A — > 00, Ui ~ 
the fh we write 



n 



2n+l 



R(X) 



-dX 



(6.2) 



with the Pij determined via / Q 



'2n+i = (again g equations in g unknowns for each 
n). Note here as A — > 00, fi 2 „+i ~ X n ~zdX so for s = In + 1 and k 2 ~ A, Q s ~ 
\ n ~2d\ ~ <iA n+ 2 ~ dk 2n+1 ~ c/A; s . Also note that the lower order terms in ( |6.2|) are 
led by (X 9 ^ 1 / X g+ ^)dX ~ A~2 ~ d(l/k) so the asymptotic behavior is correct. Then 
one determines f2 s by normalization instead of using a pseudodifferential operator 
(kl) + = k\ + 0(l/fci) to produce correct negative powers when g = 0. In particular 
p = J P Qi with dp/dP ~ f2i formally in some appropriate sense and a choice Q s ~ 
Q2n+i = J P ^2n+i will yield for s odd, dQ s /dP ~ fi s , with 



9p 



d p Qs 



gp J L s 

dp 
dP 



(6.3) 



formally (which corresponds to the result suggested in [fig, ^3| where Q s ~ &rj_). We 
note also that the branch points A^ correspond to integrals of motion (cf. ||102|| ) and 
one has Riemann invariants A& satisfying (cf. [p5 , |39| , [461 , |51 , |60| , Sections 5.1 and 
5.3, and Section 6.2 to follow) 



dXk _ f2 s . . dXk 
W s " ~?V k) dX 



(6.4) 



The Hurwitz space here is indicated in the Examples 5.1 (cf. also Section 7.3). 



6.2 Another look at KdV averaging 



we will 
1. 
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In order to better understand these matters regarding equations such as 
go over some averaging procedures from |25|] with attention to [^, |6^ 
We think first of KdV with R(X) = [i 2 = Il? m+1 (A - A fc ) and recall (cf. Section 3.1 
that for periodic situations with wave functions satisfying ( |3.2|) (i.e. Lip = Xip; L - 
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—d 2 + <fi; dtip = Aip; A = 4<9 3 — G<pd — 3<p x ) one defines the quasi-momentum and 
quasi-energy via 

p(A) = -iTtogWxi E{\) = -iJfo^Tt (6.5) 

Recall that the notation < , > x simply means ^-averaging (or ergodic averaging) and 
(logip) x ~< (logip) x > x ^ here since e.g. (logip) x is not bounded. Observe that 
( p.5p applies to any finite zone quasi-periodic situation. Now in the notation of P5" 



one takes Ai > • • • > A 2m +i with spectral bands [A 2m +i, A 2m ], • • • , [Ai, oo) on the real 
axis or m gaps (A 2m , A 2m _i), • ■ ■ , (A 2 , Ai) in the spectrum (see e.g. [§, [jjj for pictures). 
Thus (cf. (pD) 

p(A) = / dp(X) = [ V = \ m + £ (6.6) 



r r 6A m+1 -I- f(\) m 2m+i 

E(X) = / dE{\) = / ^ l } dX; £ = Y / b J \ m -i; b = -3 £ A, 

; ' ■ 'i?(A o 1 



and the normalizations are 

A2l ~ X dp(A) = f dE(X) = 0; i = 1, • • • , m (6.7) 

Note as indicated in Remark 3.4 that one will often require normalizations S J c dVL^ = 
for c G Hi(T g ,Z). Here the dj cycles can be drawn above or around the gaps 
(A 2m , A 2m _i), • • • , (A 2 , Ai) and for Uj a basis of holomorphic differentials (j = 1, • • • , m - 
g) one has (cf. (|6~ll )) 



in 



CinX^dX 



2n5 jk (j, k = 1, • • • , m); Uj = £ n / (6.8) 

R(X) 



with iBjk = § b Uj (j,k = 1, • ■ ■ , m) (the notation changes slightly from time to time in 
this paper but subsections are consistent). The matrix of periods (Bjk) is symmetric, 
real, and positive definite and the Riemann theta function is defined by 

6(t\B)= Yl exp{~J2 B jknjn k + iY,njTj) (6.9) 

— oo<ni<--<n m <oo j t k j 

and finite zone solutions <p(x, t) of the KdV equation <p t = 6<jxf) x — <fi xxx determined 
by Ai, • • - , A 2m+ i have the form 

(f>{x, t) = -2d 2 Jog9(kx + ut + t\B) + c (6.10) 
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where r — (n, • ■ • , r m ) and 

fcj = j> dp; ujj = j> dE (j — 1, • • • , m); c = ^2 Xi — 2^2 <t> \uo q (6-H) 

i 

so the branch points Aj of T m parametrize the invariant tori ( |6.10| ) of KdV (note 
kiX + LO{t + Tj ~ ^ in Remark 3.3 and note that u;.,- is bad notation here). 

Another set of parameters is provided by the Kruskal integrals Jo, • • • , li m which 
arise via a generating function 

p(A) = -i < (ZttfVO, >,= v^A + £ (2 ^- s )2s+1 (6.12) 

where J s =< P, >,= P (s = 0, 1, • ■ ■) with -i{log^) x = V\ + Eo^/^VA) 2 ^ 1 ]. 
Similarly 

- *(Zo^); = -M. = 4(VX) 3 + V § (6.13) 

V V o (2v/A) 2s+1 

and one knows <9 t P s = <9 X .P S since (Jfc) [(Zo^V)Jt = [(^°5 l ' ? /')t]a:- The expansions are 
standard (cf. || J2|, [15], |Tj|). For completeness we recall from [^| that (modulo total 
derivatives) 

P, = 0; Pi = y; P2 = f - 2 ; • ■ •; E = 3<p 2 ; (6.14) 

3 9 1 

Pi = 20 3 + -0 2 ; P 2 = -0 4 - x </w + -4> 2 XX - 34> 2 4> xx + 600 2 ; • • • 



Now consider a "weakly deformed" soliton lattice of the form ( |6.9| ) with the A» (i = 
1, • ■ ■ , 2m + 1) (or equivalently the parameters u l = If, i = 0, • • ■ , 2m) slowly varying 
functions of x, t (e.g. u l = u l (X,T), X = ex, T = et, i = 0,1,- ■■,2m). Now one 
wants to obtain a version of ( |3.17|) directly via (Jit). Thus insert the slow variables in 
(Jft) and average, using edx or e&T in the external derivatives, to obtain 



d T (logil)) x = d x (logip) t (6.15) 

or 8tp(X) = dxE(X). Then from (|6.6j ) differentiating in A one gets 

T dp = d x dE (6.16) 

which is equivalent to ( |3.17| ). Now (recall d t P s = d x E s ) expanding (|6.16j ) in powers 
of (x/A) -1 one obtains the slow modulation equations in the form 

d T u s = d x E7 s (s = 0, • ■ • , 2m) (6.17) 
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where E s is a function of the u % (0 < i < 2m). This leads to equations 

d T \ = Vi(\i,---,\ 2m+ i)d x \i (i = 1, • ■ ■ ,2m + 1) (6.18) 

for the branch points A& as Riemann invariants. The characteristic "velocities" have 
the form 

dE 6A™ +1 + £(A;) , , , 

V < = * U = 2 P(A.) (1 ^ < 2™ + 1) (6.19) 

To see this simply multiply (|6.16| ) by (A — A^) 3 / 2 and pass to limits as A — ► Aj. These 
equations (|6.18|) - ( |6.19|) correspond to (|6.4|) . 

6.3 Formulas for kernels 

We will begin by listing various formulas involving Cauchy type kernels on a Riemann 
surface. Properties and origins will be indicated along with a sketch of proofs and 
references. Then we want to examine relations between some of these kernels. One 
motivation for this investigation is the result of |T5[ which shows that the dispersion- 



less differential Fay identity is equivalent to the kernel expansion in dKP theory 

£(A, n) = p(/i) - p(A) = £ d P Q n (\)fi- n (6.20) 

(cf. Section 4). To recapitulate one can begin with a Lax operator L = d + 
u n+ id~ n , and pass to dispersionless limits via formulas ip = exp[(l/e)S(T, A)], T = 
(T n ), T n = et n ) to obtain 

00 00 

A = P + £ U n+1 p- n ; P = A - £ P i+1 A-' (6.21) 

1 1 

The KP hierarchy then becomes the dKP hierarchy in the form d n P = dB n where 
B n = A™ and P = dxS. Setting Q n = (l/n)B n one arrives at ( |6.20|) , whose impor- 



tance was first indicated in (see Section 4 and cf. also jl^, 0). In the reduction to 
NKdV the dispersionless situation involves A^ = A^ ~ E in ( |3.75|) (with P replacing 
p). In particular the action of /C(A, fi) at fi = 00 is expressed via 

1 r u n 

' ' -dfi = d P Q n+1 {X) (6.22) 



2ni Joe P(/i) - P(A) 

It was stated in JHJ that K, represents a dispersionless limit of the Fay prime form 
which we want to clarify and make more precise here. 



First we refer to [H91 pO, 51 1 where some interesting kernel results appear. To 



44 



begin with we consider mainly hyperelliptic situations for illustration and think 
°f A* 2 — lli £ ' +1 (^ — \) with branch points Ai, • ■ ■ , A 29+ i, oo. Thus, following HBT 



take Ai < ••• < A2 9 +i < oo with spectral bands [Ai, A2], • • • , [A2 9 +i, oo) and gaps 
(A2, A3), • • • , (A 2 £, A 2s+ i). We will write —A = k 2 here (cf. Q6.2p where k 2 ~ A) with 

Q s = d{k s ) + Qs m k~ m ' X dk (6.23) 

m>l 

where f2 2s = d((— X) B ) for T g hyperelliptic. Further set 

n s = (U s ) f , u 3 = £ q^k^dk (6.24) 

m>l 

(cf. (|6.1| )) and the Riemann bilinear relations imply 



Qij = Qji\ = -2iriq» (6.25) 

Also Qki = if either k of I is even and in the present situation Ss(Bij) is positive 
definite (Bij = § b u>i). Note here A belongs to a two sheeted covering of the k plane. 

• The above constructions are the same as before modulo notation but it is ap- 
proptiate now to write the quasi-momentum dp via 

dp ~ —idk; 3 f dp = S f dp = (6.26) 

Ja 3 Jbj 

This corresponds to the real Whitham normalization discussed in Remark 3.4. Thus 
Qi ~ dk and dp ~ —idk here. • 

Now we recall the (Fay or Klein) prime form (cf. |57|, |38|, |85|, 0]), which as used in 
p9| , |50| , |5T|j , has the properties: E(P, Q), P,Q G T g , is a holomorphic (— l/2)-form in 
P, Q with the properties: (G) E(P, Q) — if and only if P = Q (H) Given t a local 
coordinate on T g near Q, if P — > Q one has (cf. || |2B|, |56|, for half differentials) 

= ( f (P)- t (g))[i + o((t(P)-t(g))»)] 

y/dt(P)dt(Q) 

(I) -E(P, Q) is multivalued with periodicity conditions: ^(P+a^,, Q) = -E'(-P, Q), E(P+ 
bk, Q) = ±E(P, Q)exp(—iBkk + 2ni Jp up.) where a*,, are the basic cycles (J) Using 
local coordinates k^ 1 at 00 one takes Q s as usual with § a _ Vt s = 0. Further (shifting 
to local coordinates) 

00 

d c d z log E((, z) = -J2 tt s (()z s - l dz (6.28) 
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which implies 

00 s 

d c log E((, z) = - y £ «.(0 - + /(C) (6.29) 

1 s 

(K) In the neighborhood of oo one has (z ~ 1/k etc.) 

log^^- = E (* m + z ' m ) — + E z m z' n ^ (Q mn = Q nm ) (6.30) 

Ms = d(^) - E QsmZ^dz (s > 1); (6.31) 

Z m>l 

d3 = — U s z s ~ 1 dz; (U s ) m = / 

where the terms with Q m0 correspond to some normalization and Q ~ (ui, ■ ■ ■ ,u g ) 
represents the standard holomorphic differentials (cf. p0|). Note that the first equa- 
tion in (|6.31|) is equivalent to ( |6.23| ) since d(l/k) = — {l/k 2 )dk but in the second 
equation the (U s ) m differ from ( |6.24| ) by a factor of 2m. The Q m0 terms are some- 
times omitted in ( |6.30|) (cf. below and 0, p. 181). Equation ( |6.30|) implies now 

dz 

djog E(z, z') - —— i = E Qmoz^dz + (6.32) 

Z Z m>2 

oo in oo In i 

+ E —QmnZ m - l dz=Y J QmOZ m - 1 dz + Y,—[d{-)-^n{z)} 
m,n=l 11 m>2 1 " 41 



and consequently 



dz 00 z /n 

d z logE(z,z') = E Qmo^ m " 1 ^ + E— ( 6 - 33 ) 



m>2 Z 1 71 



Then with such terms one has /(C) = (dC/C) + E m > 2 QmoC™ -1 ^ in ( |6T29|) . We refer 
here also to (|6.74j) where it is shown that Z v {x) = d x log E(x,p) has a pole of residue 



one at x = p and f(p) = (l/2m) J au f(x)Z p (x) for / holomorphic in a neighborhood 
U of p (i.e. Z p {x) is a local Cauchy kernel). 



Now we recall that fl6.31| ) = ( |6.23| ) (up to a factor of 2m' in the (U s ) m ) and in (|6.7 ) 
one sees that Qj = j P Qj (with p = J P Q\) leads formally to dpQj = (Qj/Qi) with 



in the form djXf. = — (Qj/Qi)dx\k = —dpQjdxXk- On the other hand in the 
LG theory for genus as in Section 3.10 one has (using Qj now for Qj in Section 3.10) 



Qj ~ k + ~ <9,-S and via (|3.72| ), djE = —Exd p Qj at p = q s (E(q s , T) = E s corresponds 
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to Riemann invariants as do the branch points Xj in (|6.4j) and we shift back X — > —X 
now) Now E ~ k n so djE = nk n ~ l djk and dxE = nk n ~ 1 dxk implying also 

djk(q s ) = -d p Vtjd x k(q s ) (6.34) 

(with complex Riemann invariants k s = k(q s )) and Qj ~ jQj where Qj corresponds 
to the Qj of dKP theory in Section 4 (cf. also (|6.20|) and ( |6.22| )). We recall the time 
variable changes in Section 3.12 and in Section 4, so one can simply ignore or adjust 
the factor of j arising in dj for the moment. Then ( |6.34| ) or (|3.72|) compares with (|6.4| ) 



using the the analogy ttj/Qi ~ d p Qj ~ dpQj in a suitable time scale. This says that 
the finite zone (hyperelliptic) Riemann surface theory with averaging yields branch 
point equations ( |6.4| ) for A& whereas the LG theory for E of Section 3.10 in genus 
zero yields corresponding equations ( |3.72|) for Riemann invariants E s or k s relative to 
the n-sheeted Riemann surface of E (the number of branch points may differ of 
course). We refer to Section 5 for clarification of roles and will discuss this all later at 



some length. The formula analogous to the kernel expansion (|6.20| ) is then (cf. 
where such formulas were first suggested) 



£(A,A*) ~E?r( A )/^ or ^^)-T,7fW— (6-35) 
i "i i "i J 

depending on how we scale the Tj. Now from (|6.28|) we have d^d z log E((, z) = 

5 = -ET^jiO^^dz and we could also take d c logE((,z) E?n s (()z s /s + 

(d(/(). Note here £ ~ 1/A and z ~ with say A, \i — ► oo so from ( p.31| ), with 
abuse of notation, Q S {C) = d(l/C) - QsU^dC, ~ d(\ s ) + Ef ' QsmX'^dX = 
Q S (X) (the latter form for Q s corresponds to A — > oo). Hence via d^log E((, z) = 
d\log E(l/X, l//i), etc. we get 

1 l 00 

d^dzlog E((, z) ~ d\d^log E(—, — ) — Qj(X)fi~ : '~ 1 diJ, (6.36) 

A M , 

d c log E((, z) ~ to £(-, -) = - £ n.(A)^- - ^ (6.37) 

A jJ, i S A 



Now (|6.35|) gives two analogues of (|6.20| ) corresponding to ( |6.36|) and ( |6.371 ). Since 



d^djog E(( } z) has a second order pole at C, = z (or A = /i) we are led reject it and 
to suggest as an analogue of K the formula (A ~ 1/C, H ~ 1/z) 

1 ~ 00 Q u~ j 

P(/i)-P(A) J 
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As mentioned above in Q one writes 

l og ^<A = y C m z nQrnn (og) 

dC 00 z 11 

d c logE((,z)~-^+ £ r^-QmndC-Z^Q (6.40) 

S Z m,n=l n 

so Z Z (C) nas the appropriate singularity at z = ( (cf. also (|6.74jj ) Let us summarize 
heuristically in 

THEOREM 6.1. For A ~ (1/C), H ~ (1/z) one has an argument as in ( gjg ) 
or (|6.34D relating /C = l/(P(/i) — -P(A)) to a Riemann surface object K in ( |6.35|) 



for general Riemann surfaces (with suitable scaling of time variables), which in turn 
can be expressed via Z z (() as indicated. The Qj will be differentials normalized 
via § Qj = when g 7^ (depending on slow variables Tj), or represented by 
suitable powers of a LG polynomial for the dispersionless situation where g = 
(or in other LG situations with g = by suitable powers of a LG function). Thus 
formally (d(/() - Z z (() ~ fii(A)£(A,/x) ~ [dP(A)/(P(/x) - P(A))] and we note that 
(dC/C) -Z Z (C) ~ (dC/0 - ( rf C/(C--2)) ~ -(rfA/(A-/i)). Below we will also indicate 
a somewhat different kind of connection of /C to the kernel u based on ipip*. • 

Generally one appears to have a theory here involving Riemann surfaces alone. 
Averaging procedures have nothing to do a priori with the construction of differentials 
(but see here Theorem 7.3). Thus start with a Riemann surface S 9 and form a KP 
equation as at the beginning of Section 2 via a BA function, which defines flow 
variables x, y, t, • • • ~ t n . Then average as in Section 3.6 to get equations (ft.54Q , and 



proceed to general CIa as in Section 3.7 for say M. g ,i- No LG function is needed a 
priori (branch points are determined for Ai 9t \ via p and additional times arise via 
Ta, A ~ (at,i) for M. 9 n). Suppose first, for illustration, that E 9 is hyperelliptic as 
in Section 6.1, leading to ( |6.18|) - (|6.19|) or equivalently to (KM ). This means that 
the modulation equations determine deformations of the Riemann surface via branch 
point equations using natural flow parameters determined from the Riemann surface. 
However KP or KdV has only been used as an intermediate step and disappears at the 
end except for the slow variables. The modulation equations only involve differentials 
based on k± and slow variables (other intrinsic time variables and differentials arise 
for additional punctures). Thus (for one puncture) one creates a "generic" point 
Pi ~ 00 and builds the Qa and Ta leading to flS.4p . If one starts with a LG potential 
E as in (|3.68 ) to generate an n-sheeted Riemann surface of genus the matter seems 



somewhat contrived but in fact here, and especially in the case of general E as in 
( |3.75| ), many interesting things occur. A more detailed study of this situation ( |3.75| ) 
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is indicated, going beyond the description of Hurwitz spaces and Frobenius manifolds 
in Section 5. One desideratum here is to make more explicit the relations between 
the LG potential and the Riemann surface in the Hurwitz space approach. Other 
questions also arise. For example let S 9 be given, say hyperelliptic. If E ff arises from 
the Bloch spectrum of a KdV situation then one has developments as in Section 6.2; 
the wave function for KdV will apparently be the BA function for E s . However if we 
choose a BA function i/j for E 5 as indicated in Section 2, it will correspond in general 
to a KP situation. It seems clear that for some divisor D, ip will be a KdV wave 
function corresponding to a potential u whose Bloch spectrum gives rise to E 5 . It 
might be worthwhile to make such matters precise and explicit. 



6.4 The kernel u 

Now go back to the BA function for KP as in ( [2.1|) so that ip{t,k) is meromorphic 
in r 9 /oo with simple poles at P%,---,P g and no other singularities in r 9 /oo (here 
t = (ti) with ti = x, t 2 = y, t 3 = t, etc.). We use k~ x as the local coordinate at 
oo so if) ~ exp[£(t, k)) ■ (1 + J2T Xik" 1 ) for \k\ large where ^(t, k) = J2T t ik i - For 
t 2 i = this is a KdV situation. The BA conjugate differential ip^(t,[j,) is apparently 
defined in [49], [50], ^] to have simple zeros at Pi,---,P g , to be a holomorphic 1- 



differential in r s /oo, and to have an essential singularity at oo of the form ip> ~ 
exp[—^(t, k)](l + J2T° x\k~ l )dk (this definition is discussed later in more detail). Then 
it is "known" that 

' ^p(k 1 t)^(k,i ; )dk = (6.41) 



ic 

for C a small contour around oo (Hirota bilinear identity) and 

/oo _ 
ip(k, t)i)\k\ t)dx = 2ni5(k - k') (6.42) 
-oo 

for Qp(k) = Qp(k'). In ( |6.42| ) the contour could apparently be a closed curve through 
oo for example (which would become a straight line in the - degenerate - scatter- 
ing situation). This point is however consistently overlooked and should be further 
examined (cf. also |H| for KP-1 and the scattering situation). Next the Cauchy- 
Baker-Akhiezer (CBA) kernel oj(k, k',t) is defined via (we give an expanded form in 

am) 

i r x 

u(k, k', x,y,t, ...) = — / ip(k,x',y,t,- ■ -)^(k',x',y,t,- ■ -)dx' (6.43) 
According to [49 , |50| , 51] this kernel is to have the following properties: (A) u(k, k', i) 



is a function in k and a one form in k' (B) uu(k, is meromorphic in k in T g /oo 
with simple poles at Pi, ■ • • , P g , k' (C) As a function of k', uj is meromorphic in T/oo 
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with one pole k and zeros at P\ , ■ ■ ■ , P g 
(E) co(k,k',i) = 0{exp[-£{k ',*)]) as k' 
k — > k! . We will express ijy below as a suitable multiple of ip 



D) Lu(k,k',t) = 0(exp[f(fc,i)]) as k -> oo 
-> oo (F) u(k,k',T) ~ (dk'/2%i(k' - k)) as 



The material on the prime form is well discussed in |85|] for example (cf. 
also |57|] and Section 6.5 below) but the proofs in |49|, |5(| |51| of (|6.42| ), (|6.43| ), and 
properties (A) - (F) of the CBA kernel are somewhat unclear so we will give some 
discussion in various contexts. Formula ( |6.41 ) with ip* in place of ij)> (Hirota bilinear 
identity) can be proved in various traditional manners (cf. [f|, |10|, P2| ) so we omit 
comment. Now look at the integrands in ( |6.42|) and ( |6.43|) for large k, kl and Qp(k) = 
Qfp(fc'), namely, ip(k,i)^(k',i) ~ exp[J2t n (k n -k' n )}(l+EXik^)(l+Ex]k'- j )dk' . As 
x = ti varies one has a multiplier exp[x(k — k')} with \exp[x(k — k')}\ = exp[x^t(k — k')} 
(note dp ~ — idk so for large k, k',Qp(k') = Qp(k) ~ 3ft = 3lk'). Thus exp[x(k — 
k')] ~ exp[ix(Qk — Qk')] and fl6.42| ) has a Fourier intergal flavor. On the other 
hand from d n ip = B n i/j and d n ip^ = —B^ijj* (true at least for ^ ~ if)*) one obtains 
(Sp(fc') = $sp(k)) 



d n / ipijj^dx 



[{B n ^ - il>(B*^)]dx = 



(6.44) 



provided integration by parts is permitted. This is not obviously true but one can 
assume it under reasonable circumstances. Further for large t n and k' ^ k, \ip4>*\ ~ 
0(\exp[t n (k n — k' n )}\) = exp[t n ?R;(k n — k' n )} — > for some n when t n — > oo or — oo 
(with < isp{k') = Qp(k) or not). This (with ( |6.44|) ) implies that the integral in ( |6.42|) 
is for kl ^ k and Qp(k') = Qp(k). To show that (|6.42|) actually gives a delta 
function for Qp(k') = Qp(k) one is referred in [SIJ to |57|] for a discrete version whose 
proof is said to be extendible. This is not clear but see also |76|, [fjj which don't 
seem to explain much either. Finally it is clear that the ±oo limit of the integral 
in (|6.43p may change when k crosses the path ^sp(k) = Qp(k'). In this regard recall 
\exp[x(k' — k)]\ = exp[x{^tk' — $tk)] and ^p(k) ~ — 3ftfc for k large. Hence u in ( |6.43p 
can have a jump discontinuity across the curve C : Qp(k) = Qp(k'). To see this and 
to show that u(k, k') is nevertheless continuous for k ^ k! let / C C be a small arc 
seqment and D D I a small open set with boundary 3D. Then J D dujdkdk = J dD uodk 
by Stoke's theorem. Think of / as a small straight line segment with D shrunken 
down around / to be lines above and below with little end curves. One obtains in an 
obvious notation 

dudkdk = I (u+ — uJ)dk = (6.45) 



id Ji 

via ( |6.42|) and fl6.43|) (using the change in integration limit in ( |6.43D ). The other prop- 
erties in (A) - (E) are more or less natural for suitable ij}* as indicated below. Property 
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(F) is suggested by e.g. (l/2m) J*^ ipipUx'dk' ~ [l/2m(k - k')] ■ exp[x'(k - k')] ■ 
0(l)dx'dk' when say ^tk > ^ftk', k' —>■ k and x > with the expectation here that 
l/(k — k') will emerge from the integration (cf. also Theorem 6.3). We will suggest 
formulas for u and ip^ below which accord with (|6.41|) , ( |6.42| ), and (A) - (F), but a 
further infusion of rigor would be welcome. 



6.5 Remarks on general kernels 



We refer here to ||103|| for some general information on Cauchy type kernels on a 



Riemann surface £ s (cf. also |93j). Take a Riemann surface S 9 with a canonical 
homology basis (aj,6j) 1 < i < g, and abelian differentials u k , 1 < k < g, with 
&i = Sik- The B periods are B k j = § b Uj as usual. Let u) qs {j>) (or duj q s){p)) be a 



normalized (i.e. f a tu qs = 0) differential with simple poles at q and s with residues 
±1 respectively. One knows that w qs (p) —ui qs {r) = uj pr (q) —u pr (s) so uj qs {p) is a many 
valued function of q with a single pole at q — p. Further 

uj qs (p) = 2ixi [ uj k (6.46) 



(cf. (|6.76| )). Consequently 

f d q uj qs (p) = 0; f d q u qs {p) = 2muk (6-47) 

Choose a brance oJ qs {p) now which is a single valued function on S 9 and satisfies the 
condition duj ss (p) = 0. This branch can be defined via 

u qs {p) = J d t u) ts {p) (6.48) 

where the path of integration does not intersect ai, • • • , a g . One uses + (resp. — ) to 
destinguish objects associated with the left (resp. right) bank of the oriented curves 
afc. Thus a~l (resp. a^) denotes a k described in the + (resp. — ) sense. The surface 
S 9 is on the left of the curves a^ft^afftj" • • • a+b+a~b~ ~ dT, g in the standard cutting. 
Then from ( |6.47| ) - ( |6.48| ) one has for t e a k 



rq 

&u(p) - = -27rzw fc (p); - u~(t) = -2m I u k (6.49) 



(± over t would be slightly better notation). 

Consider now a Cauchy kernel K = {dr/ij — z)) in the complex plane, which 
has the properties: (1) As a function of z, K is a meromorphic function with a 
single pole at z = t and a single zero at z — oo. (2) As a function of r, K is a 
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meromorphic differential with two simple poles at r = z and r = oo, with residues 
±1 respectively. We denote analogues of Cauchy kernels on S g by A(p, q)dp where 
(3) A(p,q)dp = (dp/(p — q)) + regular terms as p — > q. Generally for g > a 
meromorphic function cannot have one simple pole so there is no analogue of (1). 
Hence, abandoning (1), there exist infinitely many analogues of a Cauchy kernel and 
two especially are constructed in [ |103|| . The first is a discontinuous analogue which is 
fulfilled by Cj above where Cj is discontinuous in q along the a^. Thus 

A(p, q)dp = Cj qs (p) (6.50) 

(s can be fixed arbitrarily and r, z, oo ~ p, q, s). Except for the discontinuity (1) 
and (2) are satisfied. 

The second type is a meromorphic analogue. For this one says first that a divisor 
D is minimal if ord(D) = g — 1 and r{D^ 1 ) = i(D) = where the meaning of r, i is 
visible from the Riemann-Roch theorem r(D~ 1 ) = ord(D) + i(D) — g + 1. However 
note here that this is in multiplicative notation for divisors, D = p™ 1 • ■ ■ p^ k instead 
of J2 n iPi so one knows (cf. [^]) R(l) = 1, = g, and ord(D) — ~ D — 1. 
We recall also for rii > with D > 1 one has C C L^D" 1 ) (C = £(1)) where L 
is standard notation. Thus if ord(D) = g — 1 with r{D~ v ) = then D must have 
poles. Now (given a minimal divisor D) one constructs a Cauchy kernel A(p, q)dp 
having property (3) which as a function of p is a differential with divisor a multiple 
of q~ x D (~ D — q additively) and as a function of Q the kernel is meromorphic with 
divisor a multiple of p^D^ 1 . Thus pick e.g. 

k i 

D = Y[p" J ■ n^" 1 ^ 1 (« > °> fc > °> ™j > 0, > -1) (6.51) 
1 o 

with ord(D) = Y,i n j ~ Ylo m i — i — I = g — 1. By remarks above there is at least 
one q, say q , with a positive multiplicity m + 1 > 0. Then one can write down a 
formula in determinants for A(p,q)dp (cf. [ |103|| ) with entries 



(?) = V ^f; uj£(p) = -rf^o (P) (6-52) 

This kernel will be unique via minimality of its characteristic divisor. However there 
exist other meromorphic analogues having different properties (cf. below). Exam- 



ples of various sorts are indicated in [103]. Discontinuous Cauchy type kernels are 



connected to Plemlj-Sokhotskij type formulas (with eventual application to Riemann- 
Hilbert problems possible) and meromorphic Cauchy kernel analogues are related to 
kernels of Chibrikova, Gusman-Rodin, Koppelman, Tietze, Vaccaro, and Weierstrass 
for example. 
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Some other analogues of Cauchy kernels arise from the differential 

„ w + ( dr 



2C 



r — z 



(6.53) 



where (z,w) G S g , (r, Q G £ 9 when S 9 is the hyperelliptic Riemann surface of w 2 — 
f(z) = for f(z) = T\i g+1 (z — z k ) . One thinks here of a k cycles surrounding r 2 k-i, ?"2fc) 
with (z, v7) the upper sheet. One can take as usual (dz/w), • • • , (z 9 ~ l dz / w) for the 
Uk (unnormalized). Then the discontinuous analogue of the Cauchy kernel is 



duj 



(z,w),(z ,w ) 



V, C) = duj (Z:W) (r, C) - du {z0jWo) (r, () 



where 



7- \/ a\ iVwm 



Jai ( 

£ *L 
Ja g £ 



fa 

§a,_ 



Num = 



w+C dr dr 
2C ' T-z C 
w+C dr £ dr 



*g C 
rf- 1 dr 



(6.54) 



(6.55) 



ai 2C, t-z Ja i C 
'a-g 2( t-z Ja g ( 



ai 



T<>- X dT 
C 

T3' 1 dr 

c 



One could also phrase this in terms of normalized uoj. On the other hand a meromor- 
phic analogue (elementary Weierstrass function) with characteristic divisor (zq, wo) _1 ( r i ) Ci) 
has the form 



A[(z, w), (r, C)]dr = A [(z, w), (r, Q\dr - A [(z , w ), (r, ()}dr 
where A [(z,w), (r,()}dr = (Num/ Den)(dr/2() with 



(6.56) 



Den 



1 71 
1 T n 



n 



9-1 



T; 



9-1 



(6.57) 



iVwm = 



T — Z 
W + Cl 
T\—Z 

">+C 9 



T 
7-1 

7V, 



r 9-l 

9-1 



-9-1 



fo. Ca) 
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6.6 The prime form 

We will collect first further information and some calculations involving the prime 
form (cf. H [|, H; 13> HI HI H> H> EH)- I n particular we want to confirm (|6.27| ) - 



( |6.31| ). Thus work with a lattice Lb = Z 9 + BZ 9 and a torus C 9 /Lb (we think here 
of I? with positive imaginary part as period matrix, so some factors of i may need 
adjustment relative to previous notations). The theta function with characteristic 
a, (3 is defined via 

9[a, /3] (z\B) = e i7T{n+a)B{n+a)+2ni{z+ ^ (6.58) 

where the periodicity condition 8[a,j3](z + Bn + m\B) = exp(—iimBn — 2im(z + 
(3) + 2iria ■ m)9(z\B) determines 9 up to a constant. For some arbitrary base point 
Pq define now the Abel-Jacobi map 1(D) = ^2r}ifp*u for a divisor D = J^Vi^i- 
Then f(P) = 9(z + I(P)\B) either vanishes identically or it has exactly g zeros 
Pi, ■ ■ ■ , P g . In the latter case there exists a vector A ~ Riemann divisor class, ~ —K 
in Section 2, depending only on Pq and the canonical homology basis, such that 
z + I{J2-Pi) = A. The set of z G «/(S 9 ) = Jacobian variety where 6 vanishes is a 
subset of complex codimension one called the 9 divisor. Then 9(z\B) = if and 
only if there exist g — 1 points Pi, • • ■ , P g -\ in H g such that z = A — I (Pi + • ■ ■ + 
Pg-i)- The bundle of holomorphic one forms is called K, and this is the canonical 
line bundle of degree 2g — 2; spinor bundles L satisfy L 2 = K,. We note that to any 
holomorphic line bundle L on S g , with s a meromorphic section, one associates a 
bundle L s with the equivalence class D s of divisor s. Similarly for any divisor D we 
can associate a line bundle (locally build transition functions via divisors restricted to 
open sets). The relation between spin structures and theta functions is summarized 
well in recalling that the Riemann class A satisfies 2A = /C. Hence A is the divisor 
class of a spin structure So and the divisor classes of other spin structures S a ^ are 
given by D a = A—Ba—j3, (a, f3) e (|Z/Z) 2ff to which one associates a theta function 
9{a,f3}(z\B). For such theta functions one has 9[a,P)(-z\B) = (-l) 4a ^9[a, 0\(z\B) 
leading to a classification of spin structures as even or odd depending on whether 
the corresponding theta function is even or odd. For odd characteristics a, (3 we have 
9[a, (3}(0\B) = so by Riemann's vanishing theorem there are g — 1 points P« such that 
D a ,p = [12 Pi] ([D] denotes divisor class). This suggests that S a ^ has a holomorphic 
section h a ,p(z) having g — 1 zeros for s = Pi and indeed this section can be written 
as follows. Define 

f a ^(z,w) = 9[a,/3}([ u\B) (6.59) 



Then by the vanishing theorem f a $ has single zeros for z = w, z = Pi, or w = Pi 
(obvious with a moment's thought) so for both z and w in the neighborhood of one of 
the Pi, f aj p(z,w) looks like constant- (z — w)(z — Pi)(w — Pi). Differentiating this with 
respect to w at w = z one finds that the one form g a ,/3( z ) = J2 di9[a, (3](Q\B)wi(z) 
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has g — 1 double zeros for z — Pi. Note that since the divisor of g a ,/3 is in fC we 
have indeed 2D a ^ = /C and there exists a holomorphic ^-differential h a ^ such that 
9a,p = h 2 a l3 which is the holomorphic section of S a ,p desired. Now one can define a 
holomorphic differential form of weight (— 1,0) x (— |, 0) via 

E(z, w) = - fa ( A " ,W ] r ; (a, /?) odd) (6.60) 

Let us note also that one can construct arbitrary meromorphic functions or differen- 
tials via E. Thus given that D = P\ + • ■ ■ + P n — Q\ — ■ ■ ■ — Q n is the divisor of a 
meromorphic function one can express this function via (cf. f85fl ) 



f(z) = WE(z, Pl ) 

On the other hand for example a differential of third kind with first order poles at P 
and Q with residues ±1 respectively is given by 

u(z; P, Q) = d z log^^dz (6.62) 

For the prime form (|6.27 ) we note also the following representation from [j35H , Vol. 1, 
p. 160 (E is treated as a function here). 

E(x, y) ~ E e (x, y) = 9{e + f u); 9(e) = (6.63) 



Here uj ~ (uj) and we note that 6(e) = if and only if there exist points P±, ■ ■ ■ , P g -% 

" -Pi 
Pi) 



such that e = K — Yfi 1 jpl £>. 



6.7 Remarks on theta functions 

Let us now extract some information from JJS] regarding various kernel formulas and 
theta functions. We recall the prime form as in ( |6.60| ) rewritten here as ((a, (3) ~ 5) 



h s (x)h 5 (y) 



dzi 



(6.64) 



e[Q](z)=J2e 



\mBrn r -\-rnz 1 



E 



i (m+a)B(m+a) T + (z+2m/3) (m+a)' 



(6.65) 



e 2 



\aBa T ' + (z+2ml3)a 1 



6(z + e) 



55 



Note that 9(— z) = ^exp[(l/2)mBm T — mz T ] = 9(z) via m — > —m. Here e = 
2irif3 + aB = (/3, a) ( 2 ^) and the period matrix B is symmetric with 3tB < (Siegel 
left half plane). E(x,y) will have the following properties (cf. also Section 6.4 for a 
physics approach). (A) E(x,y) = —E(y,x) (recorded earlier) (B) If A = J2i a % an d 
B — are divisors then (cf. ( |6.62| )) 

dlogfl^p^ = u B . A (x) (6.66) 

where u>b-a is a differential of third kind having poles at 6, and a, with residues ±1 
respectively. One will have also for X, Y divisors of degree n 

Y i-B 

^b-a = / u Y -x (6.67) 

X J A 

(C) If £ 9 is realized as a covering of P 1 via z : S 9 — > P 1 then 

»> - + ? r m - w - } (6 - 6s) 

where rrij = (l/2n) J a . darg[(z — z(y))/z — z(x))] and paths of integration are taken 
within S 9 cut along its homology basis in a standard manner. Note here also that 
from (|6.66|) 

exp{ / u y _ x \ = — — r 6.69 

J P E(x,q)E(y,p) 

(D) For any nonsingular / G (9) (i.e. 9(f) = and (dd / dzi)(f) ^ for some i) define 

9 89 9 8 2 9 

H f( x ) = E iB/M*): = E ^T(/M(*)wi(ar) (6.70) 

1 2 i,j=l i 3 

Then one has holomorphic Prym differentials with # — 1 double zeros 

( e(y-x-f) \ 2 _ u( ^ u r Qi 

V E(x,y) J 

It follows that for / G (9) nonsingular 

0(y-x-f)6(y-x + f) . , . . 2 , 

= E(x, y)E(y, x) = -E(x, y) 6.72 

Hf{x)H f {y) 

(E) For 5 = d x d y logE{x,y) (cf. flOl ) 

S(x, y) = d^dy Zog ^(x, y)dxdy = d x d y log9(y — x — f)dxdy (6.73) 



H f (x)H f (y)exp([ X ^L) (6.71) 

Jy Hf 
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is a well defined bilinear meromorphic differential independent of the nonsingular 
/ G (9). (F) Z p (x) = (d/dx)log E(x,p) (cf. ( |6.62D ) has a pole of residue 1 at x — p 
and Zb{x) — Z a [x) = uj^ a (x) = J^E(x,y)dy. For / holomorphic in a neighborhood 
U of p one has 

so Z p is a local Cauchy kernel. 



6.8 Expressions for uj 

First note that ( |6.31| ) is the same as ( |6.23| ) for z = 1/k and the expression ( |6.28| ) is 
a standard way of obtaining a differential with a second order pole (cf. ||). Thus 
u 2 (z, C) = d ( d z logE(z, C) = d c (E z /E) = (E z( /E) - (E Z EJE 2 ) and via E z ^ 0, E c 
at z = C the last term looks like \j{z — Q 2 (note also Qj 2 dz has zero periods - 
cf. [8"5|). From this one picks up differentials of the second kind with poles of order 
n + 1 at z = ( via ui n +i = d™~ Cj 2 {z, (n — 2, • • ■) for example and we note that 

(cf. |§) 

f b cu n+1 {z, Qdz = —D^MC) (6.75) 
where u>i = fid(. Further (cf. ( |6.62|) ) 

I u(z,P,Q) = 2m ( P uo (6.76) 

Jb JQ 

where b ~ (b 1 , ■ ■ ■ , b g ) and uj ~ (a?i, • • • , uj g ) represents the standard holomorphic 
differentials. Here the relation ( |6.76| ) follows from standard bilinear identities as does 

f uj 2 dz = 27u^- (6.77) 
Jb k d( 

and (|6.75| ) follows from (|6.77|) . 

We can make further comments now relative to ( |6.28|) - (|6.31|) by writing (z ~ 1/k) 

oo 

u> 2 (z,0 = d z d ( logE(z,() = -Y,n.{z)('- X d(i co n+1 (z,C) = (6.78) 

i 



71! 71! " \ n 



(the latter for n > 2). This leads to 

1 



£l n (z)d( = -nuj n+1 (z,0) = — — — d™ 1 uj 2 (xX)k=o (n > 1) 

[n ij. 



(6.79) 



57 



2m{U s ) m ~ i O s ~ j> oj s+1 dz\{ =0 = (6.80) 



1 r 9tt? 

-of 1 f M*lc=o = -r^r /m(Olc=o 

Si Jbm Si 



This implies in particular that (1/ s\)(U s ) m ~ Di / m (C)lc=o- We note also that 
differentiating (|6.30|) one obtains 

d c d z (log^^-) = d ( d z logE((, z) - d s (-^-) = (6.81) 
C-z (-z 

oo n—1 oo oo 

= d c d z log -y^rdzdC = E E C" V"^^ 

1 s 11 

while (fTgD and flOID imply 

OO 1 oo 

^4/0^ £ = - E ^ -1 ^M(tj) - E Q S mC m_1 rfC] = (6-82) 

1 <"> 1 

00 oo oo 

= E w'^c— 1 + E E g-m^-'C"* -1 ]^^ 

1 1 1 

One sees that these equations are consistent and the terms involving Q m0 in ( |6.3(J[ ) 
play the role of integration "constants" used for normalization in |5(J (they are absent 
in I). 

Now consider from QOSp d x tu(k,k',£) = {l/2m)ip(k,t)if) i (k , ,tj (see below for re- 
lations between ip* and ^ and further discussion relative to the definition of uj). 
Thus first let us go carefully over the construction of BA functions etc. in Section 
2. In Q2.1D for example we pick O 1 ~ 0\ = dk + • • • , J ~ Oj = d(ki) + • • • with 
§ a , flj = (or sometimes 3? § Oj = = ^§ bi Oj as in Section 3.2 for example). Set 
Uj = fa Q\, Vj = fa 2 , " " " U P to factors of 2m and recall dp ~ Oi, dE ~ 2 (or 
f2 3 for KdV). The flow variables arise via q(k) in Section 2 but the Riemann surface 
contributes via the argument xU + yV + yW in the theta function to establish linear 
flows on the Jacobian We recall a few additional facts about BA functions 

etc. from f27|, |33 . First from (12.11) we have 



^(/V ^' + fV' (6.83, 
Jpo 9\A(P) + zo) 

for z = -A(D)-K (O ~ + - • • , f7 ~ • •)• We recall also that F(P) = 

9{A{P) — e) is either identically zero or F(P) has exactly g zeros including multiplici- 
ties on the cut surface Y, g . Note here e = (ei, • • • , e g ) and A(P) = (Jp o Wi,"',Jp u) g ). 
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Then if F(P) ^ and Pi, • • - , P s are its zeros one has A(Pt, •••,P g ) = A(Pi) + 
■■■ + A{P g ) =e-Kon J(E g ) where A(P ir --,P g ) = £? (jg • • • u g ) = 
(Ei fpo • • • j Ei JrJ ^g)- Recall here that K is the vector of Riemann constants 
and 2K = — A(K^) where Ky, ~ equivalence class of meromorphic differentials. 
Hence the function 9(A(P) — e) is identically zero if and only if e can be written as 
e = Ei A(Qi) + K where Q = E? Qi is a special divisor (i.e. there exists a noncon- 
stant meromorphic function with poles only at the Qi); then for D = EPj> A(D) = 
e — K ~ A(D) = A(J2Qi) — for any D. A divisor in general position is not 
special. The above leads to (e = £ + K, ( arbitrary) A(Pi, • • • , P g ) = ( and further 
the zeros e of 9(e) = have a representation e = A(P\) + ■ • • + A(P g -i) + K where 
Pi, • ■ • , Pg-i are arbitrary (cf. [0). In particular one sees that the poles of the BA 
function (|6.83 ) arising from zeros of the denominator lie at the points Pi (i.e. for D 



nonspecial, D = E Pi, 9(A(P) — A(D) — K) has exactly g zeros Pi, • • • , P< 



For completeness let us also include here some of the formulation of (cf. 
also 0|) but omitting any conformal field theory. One writes £(t,z) = Ei* ' t n z n with 
ty(t,z) ~ exp(-0(l+Ef 'w k z~ k ) ~ exp(-£)T(t+[z\)/T(t) where t+[z\ ~ (tj+l/jz j ). 
Similarly V(t,z) ~ exp(l)(l + E?° w k z~ k ) ~ exp(£)r(t- [z])/r(t) (* is used in gjj). 
Thus ^ ~ our normal ■?/> and ^ ~ our normal ^* but poles and zeros must be ex- 
amined. Then on a Riemann surface £ 9 one takes holomorphic differentials (Abelian 
differentials of the first kind) in the form uij = ujjdz with § a .oJj = dij, § h ujj = 
Bij, Uj(z)dz = rf(E«>o^"/ n ) = — 2 P n z~ n ~ l dz (we inserted a plus sign in Uj). 
For differentials of the second kind one writes (Q G E, z(Q) = oo) 



u n Q = uo n Q dz; f u n Q = 0; / u% = 2mP n ; (6.84) 



If we write z — 1/C, this becomes cjq ~ d(l/£ n ) — E <?nmC m "C which agrees with Q n 



in ( 6.31 ) for q nm = Q nm . Further 



d z d w log E(z,w) = [- — + g„ m 2 n 1 w m x \dzdw (6.85) 

(z — w) m ,n>0 

which agrees with fl6.82| ) (since l/(^-w) 2 = (d/dz)(l/(w—z) = (d/dz) Eo° z n w~ n ~ 1 = 
Ei° nz n_1 w _n_1 ). The Szego kernel is written now as 

«,(,,„) = ^-ff^ 1 ^ =— !— + E (6.86, 
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where c G C 9 with 6(c\B) ^ (note the similarity to ( |6.40| ) which however is not 
equivalent). We will not consider Szego kernels further (nor deal with Krichever- 
Novikov kernels as in [IS, 77]). Here the Abel- Jacobi map is written as 

JQ n 1 1 

Here we have inserted a factor of 1/n in (|6.87|) , apparently forgotten in Note 
our ujj ~ z~ 2 dz whereas in ( |6.1| ) or ( |6.8j ) for hyperelliptic situations &j ~ A _3//2 <Pa ~ 
2k~ 2 dk for A = k 2 so the growth condition matches. Finally set <p°(z) = 1 with 



E— z~ m = f "3 K~*U (6.88) 



m 

<j) n (z + oi) = n (^); n (^ + 6j) = <p n (z) + 2ttz/£ 
We record now the tau function via 

T e (t, E g ) = e^6(I(t) +c\B); (6.89) 

oo 

?(0 = E QnmtJm] I(t) = (P(t) = (£ I#„) 
m,n>0 1 

and the BA functions have the form indicated above, or more generally, in the form 
required here 

* ~ /(g )e-Sr^W g ^ + _ /(g)+c|fl) ; (6.90) 
V ; 6(I{t) + c\B) V ; 

JW 0(J(t) + c|S) 



where f(z)ydz = [ydzwj ' E{z,w))\ w ^ 00 = [1 + 0(2 1 )]vcfe. Poles and zeros of ^ 
should however be examined. For comparison to ip, ip* evidently must correspond 
to U — xU + yV H — •, and I(z) ~ — A(P) (i.e. 2; ~ P and perhaps we should restore 
the minus sign here in ujj). In this spirit ( |6.83 ) corresponds to \I/ via 

and as P — > Po (z — * 00) one has ■?/> ~ \&// as required. We will not pursue this. 

Now look at u in (|6.43|) . One can take ip as in (|6.83|) and for ip* we know it is 
determined via a dual divisor D* (cf. [[TpJ, p9"| ) such that D + P* is the null divisor 
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of some meromorphic differential Q = dk + a(dk/k 2 ) + ■ ■ -. Thus D + D* — 2Q ~ 
-Ke (Q ~ oo, 2if ~ -A(Kx) so A(£>*) - A(Q) +K = ~[A(D) - A(Q) + K\. The 
prescription of ip^ in (A) - (F) seems to require for example via ( |2.2j ) (P ~ k, P' ~ A;') 

^(Jf, t) ~ e(A(P') + 2b)0(^(-PO + 4W{k\ t )dk' (6.92) 

(modulo some normalization and adjustment for path independence). We recall now 
that by definition of ip*, ip(k', t )i/j*(k', t )Q(k') = $(/c', t ) is a meromorphic differential 
in k! having its only (double) pole at oo. Given i/j and ip* as in (|2.1|) - (|2.2|) this means 
we could take Q in the form 8(A(P') + zo)0(A(P') + Zq)Q where fl is any meromorphic 
differential with its only pole at oo of the type indicated. But this is exactly what we 
achieve via ( |6.92| ) except for the pole at oo and the question of global definition. We 



recall from ( [2.1|) and remarks after Theorem 6.2 the changes induced by a path change 
7 = Ein k a k + Trirrijbj. Thus (cf. g, p. 56) one has J^Q Jp o fi+ < M,U > 
and A(P) -> A(P) + 2mN + BM with 6(z + 2ixiN + BM) = 6 (z)exp(- (1/2) < 
BM, M > — < M, z >). Thus introducing theta functions to cancel poles etc. leads 
to problems of global definition and we are better advised to use the prime form E. 
If we take E in the form E(x,y) = 9(e + uj) for example (cf. ( |6.63|) ) then one 
has, generically, for suitable paths, upon cutting the Riemann surface in a standard 
manner, 6(e + fp.uj) = 6(e + A(P) — A(Pi)) and this is unchanged with a path 7 
(which affects both A(P) and A(Pi)). However E is multivalued as indicated after 
( |6.27| ). Then instead of ( |6.92|) , one can think of 

^{k',t) =r(k',t)f[E(P',Pjf[E{P',P*)dk' (6.93) 



(recall E is treated as a function here). Thus we think of defining 00 in ( |6.43| ) as 



oj(P, P', x, y, t, ■ ■ ■) = / ip(k, x', y,t,-- -W(k\ x', y,t,-- -^(k^dx 1 (6.94) 



2ixi J± 

(up to some factors for normalization and global definition indicated below) where tp 
and ip* are given in ( |2.1|) - ( p.2|) and 



Q(k') = f[ E(P', P^ f[ E(P', P*)dk' (6.95) 
1 1 



Then u will coincide with ( |6.43| ) when ^' is given by ( |6.93| ) and properties (A) - (F) 
will hold. Thus the kernel in ( |6.43| ) can be written tentatively via 



u x (P, P', t) = ^(x', k, y,t r - k', y,t,-- -)Q(k') = (6.96) 

2m 



61 



e (^-/; o )[^+-l e{A{P) + u + Zo )e(A(n -U + 4) Of E(P', Pj)E{P\ Ft) , 

2m ' 6(A(P) + z )6(A(P') + z* ) 

Now to determine global definition we are concerned only with separate path changes 
for P or P' integrals. Thus for P the exponential term contributes a multiplier 
exp(< M, U >) and the theta functions give exp(— < M, U >) which cancels. For P' 
we get exp(— < M, U >) from the exponential and exp(< M, U >) from the theta 
functions which again cancels. Hence ( |6.96| ) is globally well defined (but multivalued) 
and we have tentatively 

THEOREM 6.2. A path independent expression for u x having the essential 
poles and singularities stipulated in (A) - (F) can be written as in 



A definition of uj via ipip* (instead of iffif)*) seems more appropriate in dealing with 
dispersionless limits (as seen in Section 6.9) but in order to have a Cauchy kernel the 
poles from iptp* should be eliminated (hence ip^). This leads to ip* as in (|6.93|) and uj x 
as in ( |6.96| ). The pole in u at k = k! will emerge from the integration as indicated 
before (cf. also Theorem 6.3) and we will have a discontinuous Cauchy kernel analogue 
in the spirit of Section 6.5. We have to rethink ( J6.41 ), ( |6.42| ), and d n ij^ = —B* n ip^ now 



but, since no x, y, t dependence has been introduced, it seems that d n ip^ = —B^ip^ 
will still hold, (|6.41|) should be OK, and (|6.42|) should be preserved up to a multiplier 



(which we don't compute here). 



6.9 Averaging for theta functions 

In terms of averaging one has (Jp Q Qi ~ p) 



d x logu x = -^ = [p(fc) _ p (k')} + d x log6(A(P) + U + z ) + (6.97) 



+d x log6(A(P')-U + z* ) 

The averaging process should now evidently revert to integrals (l/27r) 29 / • • •/ d x log9 
~ (1/2tt) 29 /■•• JJ2Ui(d/d6i)log6d 29 6i and we recall that 9(z + 2me k ) = 6(z) and 
6(z + Be k ) = exp(-(l/2)B kk - z k )d{z) (cf. |7[). This would imply for the fourth 
term in (|6.97|) (note (9i,8 g+ i) ~ Zi), running 9j between and 2n, 



j "... J d x log6 = -f^Uii-^Bu) (6.9? 



(here 9j ~ A(P) + Uj + (zo)j and running this in (0, 2n) eliminates the zj terms 
from —(l/2)Bjj — Zj above). But from the fifth term in (|6.97|) we would obtain 
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correspondingly — J2l Ui{— (l/2)Bu) yielding formally 

(4) + (5) = -£tM-Ifl«) +j^ Ul {- l -B ll )] = (6.99) 

Note here Ui ^ Ui since U = xU + yV + • • • and Ui = xUi + yVi + • • • (althought 
this is not used here). The quantities Ba, U are then thought of as functions of 
slow variables X,Y,T, ■ ■ ■. This is perhaps somewhat cavalier but starting a period 
at some Z{ ^ will leave terms in ( |6.99| ) depending on x, y,t, • • •; recall also that (up 
to possible factors of 2ni) 

Ui=<f fii; U(=<f n[; B u = I W< ; B' u = I (6.100) 

J hi J hi J hi J hi 

Consequently one obtains 

< 8 X logu x >= p{k) - p{k') (6.101) 
which agrees with (|6.5|) and ( |6.12|) for example in the hyperelliptic case where k ~ 



6.10 Connections to differential Fay identity 

We recall from that the differential Fay identity leads to 

nUW,n) = J-a{ ^gl } = (6.io2) 

11 — X T{t) 

1 5[c TM*-») T(t + M) 3 



fjL — X r(t 

Given the equivalence of the dispersionless differential Fay identity with the kernel 
expansion ( |6.20|) (cf. [15|]) one expects fl6.102j ) to be at least formally useful in studying 
lj in ( |6.43|) since (|6.102|) implies (modulo adjustment for ip^) 

u ~ ^-d^[r(x',t n ,X)^(x',t n ,fi)] ~ -L 1 *M,j*M*) 

Z7TZ H — A T\t J 

modulo integration or normalization factors (note ( |6.103| ) implies uj ~ l/27ri(/x — A) 
as /i — > A). Hence one should be able to determine easily a dispersionless limit for a;. 
Thus as in J15], [57 express r via r ~ exp(F(T)/e 2 ) and write 

' !E*-^l£iit« (6 . 104) 

2irifi-X r (T) 
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Take logarithms to obtain then 

loguo, ~ log— - logU - A) + - V TAfi 1 - X 1 ) + (6.105) 
2ixi e 

+J 2 {F{T + e[X- 1 ]-e[^])-F{T)} 
The last term can be written as (% n ~ elementary Schur functions) 

IjeE^ -e-S^^F-F} = (6.106) 

1 oo oo 

= -IE Xn(e<9)A-" ■ E Xm (- e 9)/i- m F - F} = 

6 

-1 oo oo oo oo 

= -{E( )» + E( )™ + EE( )»( U 

e 1 1 11 



Now we have (cf. (cf. 



1 oo oo oo oo rp 

4EE( )»( ) m — EE^-y™ (6.io7) 



11 1 1 ™ 



and we recall here from IT 



£ £ ^a->- = -Mi - f ) - E = ( 6 - 108 ) 

ii nm A ^ A n 

= -log(X + logX - £ Qn(^)X~ n = log[ P{X) ~ P{ ^ \ 

1 A — fl 

(cf. also Section 4). Now try an expression u e = f(X,n)exp(R/e) which will entrain 
- + logf~ log(^) - log(ji - A) + - E W - A*) + (6.109) 

1 oo oo oo oo 

+ -{E( )»+E( )m"EE( )»( )m} 

e 1 1 11 

Multiply by e and let e — > to obtain 

00 rp oo rp 

R ~ YTAy} -\ l ) + Y —X- n - Y —n~ m (6.110) 

1 n i m 
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which leads to 



log f ~ log(^-) - logQi - A) - %( P(//) f (A) ) (6.111) 

2 717 /x — A 



This implies via jT5| and Section 4 



d x P = (ji - A) + £ -^(A~* - = P<» - P(A) (6.112) 
1 ra 

and 

-Zog[P(//) - P{\)]+log{^-) (6.113) 

27TZ 

Thus w e ~ (l/27ri)[l/(P(/i) - P(A))]exp{l/e)[5(//) - 5(A)]} and djot/u^ ~ P(/i) - 
P(A) (c^. = ec?x). Note dxS = P in the notation of JT5J implies 

d x R = d x [3(ji) - 5(A)] P ~ 3(ji) - 5(A) (6. 114) 

One can therefore state 

THEOREM 6.3. A dispersionless kernel analogous to uj can be modeled on 
( |6.103j ) to be extracted from 



(recall also if) ~ exp(S/e) and ^* ~ exp(— 1/e) in the dispersionless theory). We see 
here how /C ~ (l/[P(p) — P(A)]) arises in looking for a dispersionless form of u in 
the genus zero situation. • 

We note that 2mu x ~ ipip* corresponds to 

d x (logW x ) e ~ d x - e [S(n) - 5(A)] ~ P(/i) - P(A) (6.116) 

so results like ( |6.5| ) seem natural - i.e. the dispersionless limit agrees with averaging 
upon identification of the differently defined quantities p(A). We do not bother here 
with the factor of i which arises from a particular normalization (our A,/i above 
correspond to k', k). 

7 FURTHER PERSPECTIVES AND CONCLU- 
SIONS 



We begin with some remarks related to results in [ITR ITS 
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7.1 Connections to inverse scattering 

Now referring to KdV situations and we recall the classical KdV picture a la |9|, |12| 



15| , |2l| , p. 02|1 . Thus for ut = u'" — Quu' with ip" — uip = —k 2 ip and ip± ~ exp(±ikx) 



as x — > ±00 one has T(k)ip_(k,x) = R(k)ip + (k, x) +ip + (—k,x) where R, T represent 
reflection and transmission coefficients respectively. For ip_ = exp(—ikx + <f)(k,x)) 
this leads to 

v = -dlog(V>_) =ik- 0'; 0" - 2ik<f>' + 0' 2 = u (7.1) 
✓ - I ^ - * + 1 ^ - MT) = I ^ £ M* (7.2) 

involving n = —v n , where these n are a priori unrelated to the ra of ( |6.14| ). In the 
dKdV context with ip + ~ ip = exp(S/e) and ^_ ~ 1/;* — exp(—S/e) one obtains (cf. 
0) 

1 p 00 p 

1 - 1 - 1 V rn+1 (Ti) 
T ik 

(we have written P n ~ P n (X, • • •) here to distinguish these coefficients a priori from 
the P n in (??)). One also obtains P — ^ = ik — <$>' = v; P 2 — U = —k 2 where 
■0_ ~ exp(—S/e), <p'(x,t) ~ 4>'(X, T) + 0(e), P = dxS, etc. These relations are quite 
interesting as pointed out in Hl5| . The scattering data T determined by asymptotics, 
has meaning in some sense related to Whitham averaging on a degenerate Riemann 
surface, or to dispersionless limits, and e.g. the classical action variable log\T\ can be 
thought of as depending on the slow variables X, etc. Further, formally, from ( |7.2j ) 
and Q7.3Q we get 

P = i k -t~ ik -±J^- = (7.4) 
which indicates how the classical n are related to slow variables via P„,+i Also from 

o 

Because of the first order term in ( |7.2|) for this to be correct we must have /f^ 0!<ir = 
in ([T2"D, which corresponds to if = or J^udx = (cf. jnj]). With this 
assumption we can expand the logarithm in ( |7.5| ) to obtain via log(l — x) = x — 
(l/2)x 2 + (l/3)x 3 + --- 

00 p 1 / 00 p \ * 

-M^E^-lfclSr + ■ — (7.6) 
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(ik) 2 (ik) 3 2 \(ik) 4 (ikf 

Hence from ( |7.2| ) one obtains 



oo J — oo 



= / fodx ~ P 2 ; / fodx ~ / P x dx ~ P 3 ; (7.7) 

J —oo 

rx 1 



0= / fadx ~ P 4 - ^-P 2 ; 



OO 



2 



This gives connections of the universal coordinates P s from dKdV with classical scat- 
tering data. In [12] one uses u t — uu x + u xxx = with ^ + ~ exp(ikx + <p(k,x)) to 



obtain (modulo ± signs) 

^ ~ E /Sb ^ ~ 2" • 18 / 2n+1 rfx; (7.8) 

l^rCJ J —oo 

E SHn = d 5Hn^_ E = Q3+ 1 + ^ 1 /- ^ ^ J»y_ <^ 

OM OM 3 2 J-oo J-oo 2 

We note a few obvious comparisons. First one is dealing with Lip = Xij) in Remark 
6.2, where L = — d 2 + u say, so that A ~ k 2 with y/X ~ fc. Hence (|6.12|) corresponds 
to -i(logip) x ~ fc + Er(^/( 2fc ) 2s+1 - Compare this to (jTTJ) for example to get (for 

~ 0_) 

oo roo / oo roo / j 

- ip^-)]. = * + ^ <p~ e ^ = -* E(-i) m j -tr + r (7-9) 

since 2m cfe = for KdV (cf. 0). Hence one thinks of ~ E^Ho P s dx/(2A;) 2s+1 ) 
yielding (l/2 2s+1 ) / P s cfe ~ -i(-l) s / 2s +idx or 

P s ~ -i(-l) s 2 2s+1 2s+ i = « s 2s+1 ^H s ~ ( P s dx ~a s ( 2s+1 (7.10) 



(the latter from |T2], |2T|). We also see from ( [7.2| ) that in some sense —log(T) = 
-iJ2™(-mfZo<p2s+idx/k 2s+1 ) so that 

which it is tempting to compare with Q6.12 ) for y/X = k, i.e. p(X) ~ p(k) ~ k + 
X)o°(^/(2A;) 2s+1 ). Such a comparison would involve 

/oo 
P s dx (7.12) 
-oo 
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This suggests an heuristic theorem of the form 



THEOREM 7.1. 

corresponds to p(X) ~ 
PW ~ /» fil ^ in 



Inverse scattering data — log(T) for classical KdV on R 
in the finite zone situation where the quasi-momentum 
(with general Q\ ~ fi 1 as in (|6.23 ), (|6.79| ), etc.). Relations 



to averaging of square eigenf unctions arise via formulas of the form (|3.15|) (where 
a different normalization is used) or simply via ( |3.13| ) and (|7.10 ), i.e. one relates 
the Tj, Xj arising from square eigenfunction expansions to the H s and the H s can 
be defined via inverse scattering or via differentials. Further, a direct connection of 
classical scattering ideas to slow variables is expressed via ( |7.3|) and n ~ P n +i m 
(|7.4|) when <p'(x) ~ <j>'{X) + 0(e). Other formulas hold as indicated above. 



7.2 Relations between differentials and averaged quantities 

We recall that averaging in Section 3.5 involved 7 = ifnp* directly, whereas in Section 
3.6 expansions were used involving terms like < ■0*L 1 -0 >, < i\)*A x i\) >, etc. In 
the dispersionless limit situation one could also conceivably adapt Section 6.9 to 
the latter developments but this is not immediate. Let us write from Section 3.6 
(T ~ px + Ey + nt + a ■ ( and S = xU + yV + tW + () 

< ip*L l ip >=< ip*(-2d)ilj >= -2 < e- r (f)*[pe r <j) + Ue T (f)e} >= 

- 2[p < <f)*(p > +U < >] (7.13) 
Similarly from A 1 = —3d 2 + (3/2)u we have 

< f >= -3 < ^*[p 2 e T + 2pUe £ <p e + U 2 e T <p ee ] > + (7.14) 

3 3 

+ - < i)*u i) >= -3[p 2 < 0*0 > +2pU < <p*<p e > +U 2 < 0*0 ee >] + - < 0*u o > 

We note also a natural comparison with (13.161 ) for the formula 



00 



1 A 

We remark here that such a series is natural from asymptotic expansions but when 
if}, if)* are written in terms of theta functions it requires expansion of the theta func- 
tions in 1/A (such expansions are documented in WJ\, p. 49 for example). It is now 
natural to ask whether one can express dp, dE, dfl from (|3.52|) in more detail. Thus 
note first from (3.52) 

dE < 0*0 >= 2dp[p < 0*0 >+U < 0*0e >]; dtt < 0*0 >= (7.16) 
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3dp[p 2 < <j)*(f) > +2pU < 4>*4>0 > +U 2 < <j)*(f) e e >} - -dp < 4>*u (j) > 



which says e.g. 



dE = 2pdp + 2U< ± ^ e> dp; dQ 



< 0*0 > 



(7.17) 



3p 2 dp + 



6pU- 



< 



> 



< 



> 



3U' 



< 



> 



Given p~A + -- -,E~A 2 + - 



q 



< 



A 3 + 



> 



3 < 0*w o > 
2 <6*6> 



dp 



this is reasonable. One would 



also like to compare (|3.24j ) and the equation drdp + dxdQ + dydE = (recall 
t = (1/4)9 2 7 + 39- 1 (m 57) and Q = (3/A)d y d~ 1 -f). One expects dp = d\+-~, dE = 
2\d\ + • • • , and dQ = 3X 2 d\ + • • • so a formula like ( |7.15 ) has an order of magnitude 
compatibility with a possible identification dp ~< 7 > dX. Note here the analogy 
with (|3TT6|) etc. where < ^ >~< T > with Q x ~ - < T > (dfx/2 y /Ji) (the l/^fji 
weight factor arises from KdV notation here, linked to a hyperelliptic Riemann sur- 
face); recall also < X > is similarly connected to f2 2 . Let us assume the "ansatz" (^>): 
For some choice of standard homology basis and local coordinate dp ~< ipip* > dX. 
We have seen that a variation on this is valid for KdV situations and via changes in 
homology and local coordinate is reasonable in general. In fact it generally true via 

LEMMA 7.2. The ansatz (<» is valid for KP situations. 

Proof. We refer here to fl l] for background and notation (cf. also |16|] for dis- 
persionless genus zero situations). We write d = L + J2T cr l-^~' ? which implies -$■ = 



A + EfVjA^ 

X+J2T and {log^) x 



(this is motivated by the proof of Lemma 2.1 in 

X >A 



p = A+zr < * 



). Then {logip) x 
But s n+ i = -na\-Y^T l djo\ 



and one obtains < s n+i >= — n < o 
dp = dX- j < o) > X- j - x dX = 
^* = E~ SnX~ n . QED 



> 



En 
1 



n-1 



> 



-n < oi >. 



Thus 



dX + Ei° < Sj+i > A 3 ~ l dX =< ipip* > since 



Consequently from ( |3.52| ) we obtain 



dQ 



< i)*A l i) >dX; dE~ - < ^>*L x i> > dX 



(7.18) 



This would require < ip*L l ip >= 0(2A) which by (|7.13 ) is correct while < ^A 1 ^ > 
should be 0(3A 2 ) which is also correct. Therefore one has 



THEOREM 7.3. With the notations as above dp =< i/ji/j* > dX, dE 
ij)*L x %j) >, and dQ ~ - < i\)*A x ^ > dX. 



< 



It follows then from ( 3.24Q (holding Y or T constant) that 
COROLLARY 7.4. From Theorem 7.3 results < £ > dX 



dQ and < Q > 
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dX ~ dE; consequently - < ■0*L 1 '0 >~< Q > with - < 4>* A l ip >~< t >. 
In summary we can also state 

THEOREM 7.5. The quantity ipip* is seen to determine the Whitham hierarchy, 
the differentials dE, dfl, etc., and via u> ~ /C as in Theorem 6.4 the theory is connected 
to the dispersionless Whitham theory. 

REMARK 7.6. Since D + D* — 2oo ~ ipip* is determined by a section of 
K-z (global point of view) but it is relations based on < ip*L l ip >, < ^* A 1 ^ >, etc. 
(based on the Krichever averaging process) which reveal the "guts" of ipip* needed 
for averaging. 

7.3 The moduli space M gn 

One can formulate the LG theory modeled on A n _i (or A n ) in a genus g situation 
as follows (cf. m, 0). One considers (*) d a L = [L,L+ /n ], a ^ kn, for L = d n + 
g n _2<9 n ~ 2 + ■ • ■ + qo (nKdV). Then one averages (*) over the family of g gap solutions, 
with averaged Poisson brackets, Hamiltonian structures, etc. The space M = M gn ~ 
{parameters of g gap solutions} is the moduli space of Riemann surfaces S 9 with 
genus g, a marked point Q ~ oo, and X(P) a function with a pole at oo of order n (note 
this is not the same as M. gn of Section 3.7). The dimension of M gn is iV = 2g + n — 1 
as in Section 5 (2g + n + ■ • ■ + n m + 2m = 2g + n — 1, n = ^ — 1, m = 0). One 
takes v? = X(Pj), j = 1,---,N where dX\p. = 0. Evidently the number of critical 
points is greater than n — 1 in general so A cannot simply be a polynomial of degree 
n as in the A n _i situation. The exact form of A can vary (note e.g. even a function 
like p n + q n -2p n ~ 2 + ■ ■ ■ + qo + (cn/p) + • ■ ■ + (j3/p s ) could generate more than n — 1 
critical points). In fact A will have various natural forms illustrated by the case M\\ 
where X(p) = V(2up,uj,uj') where V is the Weierstrass function (cf. f3lj). One takes 
now a covering M gn of M gn by fixing a homology basis (a^, hi) and choosing a local 
coordinate k~ x near oo such that k n = X. Then one fixes an abelian differential dp 
such that dp = d(k + 0(1)) as A — > oo with § a . dp = 0. The flat Egoroff metric on M 
in local coordinates w- 7 has the form 



N 



(dp) 2 



ds 2 = ^2g ii (u)(du t ) 2 ] gu(u) 



—Resp 



n 



dX 



(7.19) 



i 



and the flat coordinates T 1 , • • • , T N for ds 2 are 



T 



i a 



—nRes 



X » 



dp (a = 1 



• -,n- 1); 



(7.20) 



n — a 
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T n - 1+a = ^-i pdX (a = l,---,g); T 9+n ~ 1+a = I dp (a = 1, ■ • • ,g) 

ZlTl Ja a Jb a 

Note that this is a global definition of the flat coordinates so M is unramified over 
some domain in C N . Further the metric ds 2 in the flat coordinates has the form 

< dT a , dT 13 >= S a+P ' n (1 < a, p < n - 1); (7.21) 

< dT n ~ 1+a ,dT 9+n - 1+ P >= 5 a 'P (l<a,(3<g) 
and is otherwise zero. 

Let us continue this example following [|3^, [3~3fl . Thus p — Jp dp is a multivalued 
function on S 9 with Pq chosen so that \(Pq) = in some domain on M gn . Then the 
multivalued differential pdX has the form (as k ~ A 1 /™ — > oo) 

pdX = kdX - (j2 T a k a - X + OikT 1 )^ dk; (7.22) 
I pdX = 2mT n+s - 1 ; A as {pdX) = 0; A bs {pdX) = T 9+n - l+s dX {s = l,---,g) 

Ja s 

where A as (/(P)) = f(P + a s )-f(P), A 6 .(/(P)) = f(P + b s )-f(P) for any function 
or differential. The primary differentials (perhaps multivalued) are defined via 

K = d a (pdX)\ x=c = -d a (Xdp)\ p=£ (a = 1, • • ■ ,N) (7.23) 

so that A = X(p) is the LG potential of the model. Explicitly (for k — > oo) 

a = dp a = (-k a + 0(k~ 2 )dk; I dp a = (a = 1, • • • , </) (7.24) 

are normalized differentials of the second kind with <fii = dp 1 = —dp. Next 

0n-i+a = u a (a = 1, • • • ,g) (7.25) 
are the normalized holomorphic differentials (§ u a = 2ni5 a ^) and 

4>g+n-l+ a = ^a{<y = l,---,g) (7.26) 

are holomorphic (modulo dX) multivalued differential with increments 

A ag a a = 0; A h = 5 a>p dX (7.27) 



Then the primary part of the associated hierarchy for the solution where T a ' = T c 
and dxui = 1 (j = 1, c • • • , iV) has the FFM form (from 



d Ta dp = -d x <p a (7.28) 
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or equivalently in the diagonal variables u 1 ,- ■ ■ ,u 



N 



W = -(g) |pA«* = -(;§) I* U = 1,~;N) (7.29) 



(we omit discussion of the theory involving T a '° - cf. ||32|| ). Equations ( 7.28|) for 



a = 1, ■ ■ ■ ,n — 1 correspond to the averaged Gelfand-Dickey hierarchy over g gap 
solutions. The extension for a = n, • • • , iV is necessary to construct a closed primary 
operator algebra. One has 



Nil N 



Va,p = H^es Pi -^-; c al3l (T) = -J^Resp, (7.30) 

and this primary operator algebra can be represented via relations 

<Pa<P/3 — cls'P-ydp {modulo d\) (7-31) 
This is the M gn model. There are evident connections of this approach to develop- 



ments in 0, |75[ , some of which we indicated earlier. The relation to Hurwitz spaces 



as in Section 5 is clear and one can easily fill in the correspondences. 



7.4 Whitham and Toda 



We would be remiss not to mention some recent work ( [p6[ [47], |55], |80], [83], 0, |88|, £55] 
for example) on Whitham and Toda related to iV = 2 supersymmetric gauge theories 



We sketch the development in [88| dealing with N = 2 supersymmetric SU(N) Yang- 
Mills theory but will not discuss the physics here and will omit proofs. Of concern is 
the differential 



dS 



x 



dP{x) 
dx 



-dx 



y 



on the family of hyperelliptic curves 



y 



P(x) 



A 



2N. 



P(x) 



x N + 



N-2 




(7.32) 



(7.33) 



Here u = (112, ■ • ■ , Wjv) are parameters for flat moduli and A is fixed. The spectrum 
of excitations in the theory is measured by the units 



cf dS; a Di = £ dS (1 <i < N-l) 



(7.34) 
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where (a:j,/3j), I < i < N — 1, is a standard homology basis and the on circles 
counterclockwise the cut between two neighboring branch points in the x plane. Let 
Poo, Poo be the two points at oo (i.e. xQooo) = x(poo) = oo). One can write 

N-l 

dS = dX^i + ^2 a i u i (7.35) 
i 

where u>i (I < i < N — 1) correspond to standard normalized holomorphic differen- 
tials {§ a .ojj = 5ij). Here dXoo^ is a normalized meromorphic differential of second 
kind with second order poles at p^ and p^ and § dX^i = 0. 

Now define h = y + P(x) and h = —y + P(x) and consider infinitesimal deforma- 
tions of the moduli parameters u with h or h being fixed so e.g. (d / 'duN-k)dS\h= c = 
— [x k jy)dx. Changing the u — > a = (a 1; ■ • • , aAr_i) this becomes (d/dai)dS\h =c = uji 
which implies 

where (d/dcii) fixes aj (j ^ i) and h. These correspond to a subset of classical 
Whitham equations 

d d d d d d 

~ — = — — w,-; — — uii = — — Qa\ t; — r = — — riyi (7.37) 
3a* 3 daj ' <9T A 9a, ' dT A dT B y ! 

Then one specifies the Qa to be meromorphic differentials fioo.n and fioo,n of second 
kind (n > 1) where Qoo,n (resp. fioo.n) has a pole of order n + 1 at Poo (resp. p^) 
and is holomorphic elsewhere. Let also f2oo,o (= ^oo,o) be a differential of third kind 
with simple poles at p^ and poo of residues ±1 respectively. All differentials are 
normalized via § a . fioo,n = § a - &oo,n = (n > 0) and they are determined by h (h) as 
follows. Define local coordinates z^ (zoo) via ^ — h~ l (z^ = h^ 1 ) and near for 
example 

^oo,n = t-^oo" -1 " E QmnZZ'^dZoo (n > 1); (7.38) 

m>l 

^oo,n = [^n^^oo 1 ~ E W^ST^do] > 0) 
m>l 

with analogous expressions near p^. The integrability condition (|7.37|) implies there 
exists dS such that 

d d d ~ d 

— — dS = uJi] ——dS = floo, n ] -7^-dS = floo n ] ——dS = Ogo o (7.39) 

where 1 < i < N — 1 and n > 1. One then constructs a function F(a,T,T) from dS 
via 

^ = — idS{= ^M); (7.40) 
aa, 27ri Jft 27ri 
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dF 
dT n ~ 
dF 



dF 

-ReSp^z^ dS; 7^- 



-Resp^z^ dS; 



dT 



-Res Poo log z^dS + Resp^log z^dS 



Consistency of this definition follows from (|7.37|) and the Riemann bilinear relations. 
The local behavior of dS can be described by F as 



dS 



E " W 1 + T ^ ~ E 



n>l 



n> 



dF 



.71— 1 



dzoo (near p^); 



(7.41) 



dS 



- E rc^oo 1 - T ozJ ~ E 



n>l 



8f r °° 



n>l ^- Ln 



eteoo (near p^) 



An interesting class of solutions of the Whitham hierarchy are solutions such that 
Ei' 1 a i (dF/da i ) + Y,n>vT n (dF/dT n ) + Y,n>i T n (dF/dT n ) = 2F in which case the dS 
of ( |7.39|) satisfies 

N-l 

dS = a i U i + E T « fi oo,n + E T n ^oo,n (7.42) 

1 n>0 n>l 

which is consistent with ( |7.41| ). Therefore one can reproduce the dS of ( |7.32| ) by 
setting T\ = —T\ = 1 with the other Ta = 0. The pre-potential T of iV = 2 
supersymmetric SU(N) Yang-Mills theory is then given by T = 2ixiF . For this F 
one can write after some calculation (u { = —J2m>i< 7 im z ^ ) ~ 1 dz QC> near p^, with a im 
defined analogously near p^) 



N-l 



N-l 



— E T ^ a i a i + E a i 



k>l 



+ ^ E fe^fcT 1 / + E T klTkTl + ]- E ?M T fc T * 
^ fc,Z>l fc,/>l ^ kl>l 



(7.43) 



+ 9 r o,o^o + T o E( r M T ^ + r 0,fc7fc) 
Z /fc>0 



ft 2-Kl Jp, 



— ; / ^oo,o; 

7U JBi 



(7.44) 



^o,o = -Respjog + Res Poo log Zootlocfi 

To connect this with Toda let t = (ii, • • •), i = (t\, ■ ■ ■), and n be Toda lattice 
times where n G Z. One knows from [78| that the Toda lattice hierarchy has a 
quasi periodic solution corresponding to a hyperelliptic curve as above, giving an N 
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periodic solution of the Toda chain. The moduli parameters are invariants (integrals 
of motion) of such solutions. One introduces new time variables 9 = (9±, ■ ■ ■ , 0jv-i) 
into the associated BA function as 



9(z(p) - z(D) + K - ntzpOoo) - Zipoo)) + En>l tn<?n + En>l + ^ £f ~* OjTj) 

9(z(p OQ - z(D) + K - n(z(p 00 ) - z(p 00 )) + E„>i t n a n + En>i + \ Ef _1 ^Tj) 



^( Poo ) - 2(D) + K) 



(7.45) 



0(z(p) - z(D) + K) 

Here 2 is the Abel map (~ A), D is a positive divisor of degree iV — 1, .K" is the 
Riemann constant, and <7j, Oi {% > 1) and Ti (1 < i < N — 1) are iV — 1 dimensional 
vectors (e.g. <7j = T (crj!, • • • , Ct,jv~i))- As 6> — > 0, ^ — > BA function of the ordinary 
Toda lattice. A tau function can be now introduced near p^ via 



and near p^ by 



r(t,t,n,9) 



~n V t>1 ^» r(t,t -&<,], 71+1,! 



(7.46) 



r(t, £, n, 6 1 ) 



(we recall that e.g. [z] ~ (zi, (l/2);2 2 , (l/3)z 3 , ■ ■ ■)). By matching this with ( |7.45| ) one 
obtains 

r{t,t,n,9) = e p{t > i > n > e) ■ (7.47) 

1 N-1 

0(-(n - l)^(poo) + nz(p 00 ) - z(D) + K + ^ tn&n + X ^™ + ^ X! ^ Ti ) 



n>l 



where F is a polyomial given by 



1 1 

F(t, t,n,9) = - ^ qntkU + r kihU + Qklhti 

Z k,l>l k,l>\ k,l>l 



(7.48) 



1 V , n(n-l) _ . ^ 
— X, + 2 r ° ! ° ^ * 



AT-l 



k>l 



~ n X r k,oh - (n - l)r ,fctk + X 4*fc + X 44 + d n 

fc>l k>l k>\ 
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Then following H one introduces slow time variables 

and a,i = ieOi. Then from the asymptotics of (|7.45|) and ( |7.47|) one obtains 



dlogVfa -,- r ^r) = e- 1 ]T e»dSW(p, T, f , T , a) (7.49) 

6 6 6 l€ n>0 

Zo^, I, ±) = 6" 2 £ ^("'(T, T, T , a) 

The leading order terms dS^ and are the same as dS and F in (|7.43| ) and 
one thinks now of the moduli parameters u as functions of the slow variables (ut = 
Uk{T, T, T , a). One obtains a system of modulation equations for the moduli param- 
eters and these represent the Whitham hierarchy ( |7.37| ). 
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